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1.0  Introduction 


It  has  long  been  recognized  that  in  modeling  Electro-Optical 
sensor  performance  it  is  very  difficult  to  account  for  the 
ability  of  a  human  observer  to  perform  target  detection  and 
recognition  tasks.  In  general,  our  understanding  of  the  hu¬ 
man  perception  process  is  simply  not  sufficient  to  allow  for 
reliable  modeling.  These  facts  pose  a  significant  problem 
for  the  vehicle  and  vehicle  countermeasure  designers  at 
TACOM.  Other  than  the  designers'  own  intuition,  little  guid¬ 
ance  is  available  concerning  the  relative  detectability  of  a 
new  vehicle  or  vehicle  countermeasure  design.  This  means 
that  cue  features  which  significantly  influence  the  detect¬ 
ability  of  a  new  design  may  easily  pass  unnoticed  until 
prototype  vehicles  are  actually  constructed  and  evaluated. 

The  Imaging  Sensor  Simulation  Model  (ISSM)  described  in  this 
manual  represents  a  step  toward  alleviating  the  difficulties 
associated  with  predicting  the  detectability  of  a  concept 
vehicle  design.  The  model  allows  the  image  of  a  vehicle 
positioned  in  an  arbitrary  background  to  be  displayed  as  it 
would  be  seen  if  viewed  through  a  simulated  sensor  under 
simulated  atmospheric  conditions.  Currently,  thermal  imag¬ 
ing  sensors  may  be  simulated;  the  atmospheric  models  include 
provisions  to  treat  virtually  any  natural  or  battlefield 
obscuration  condition.  Given  a  single  target/background 
input  image,  the  model  allows  the  user  to  vary  the  observa¬ 
tion  sensor,  viewing  range,  atmospheric  condition,  and 
battlefield  scenario.  It  is  hoped  that  by  allowing  concept 
vehicle  designs  to  be  viewed  in  a  simulated  battlefield 
environment,  the  designer  will  be  significantly  aided  in 
attempts  to  create  less  detectable  combat  vehicles. 

The  computer  model  itself  is  written  in  FORTRAN  IV  so  that 
it  may  be  transported  between  a  wide  variety  of  computer  in¬ 
stallations.  The  model  is  menu  driven  allowing  for  conve¬ 
nient  operation.  Because  of  the  model's  large  scope,  several 
hundred  inputs  will  commonly  be  required  to  specify  a  single 
unique  calculation.  To  help  manage  these  inputs,  an  input 
library  scheme  has  been  implemented  so  that  a  number  of 
default  sensor,  battlefield,  and  natural  atmosphere  module 
input  sets  can  be  made  available  to  the  user.  Selected  input 
sets  may  then  be  modified  from  appropriate  user  menus. 

The  model  is  currently  running  on  an  IBM-compatible  Amdahl 
Model  5860  computer  system  located  at  the  University  of 
Michigan,  and  a  Digital  Equipment  Corporation  VAX  11/750 
located  at  TACOM.  The  two  versions  of  the  model  are  func¬ 
tionally  identical;  however,  this  report  specifically  treats 
the  TACOM  installation. 


2.0  Technical  Manual 

The  overall  modeling  approach  is  depicted  in  Figure  2-1.  In 
its  principal  mode  of  operation  a  two-dimensional  scene 
radiance  map  input  to  the  model  is  transformed  into  an  image 
representative  of  a  modeled  sensor's  output  just  prior  to 
its  display.  The  units  of  the  input  radiance  map  are  W/cm^- 
sr-micron;  it  is  assumed  that  the  scene  radiance  is  constant 
over  the  spectral  bandpass  of  interest.  The  output  map  of 
sensor  signal  and  noise  has  units  of  watts?  a  voltage  map 
may  be  obtained  by  multiplying  the  output  by  the  modeled 
sensor  responsivity  (V/W)  should  this  be  desired. 

As  illustrated  in  Figure  2-1,  the  input  scene  radiance  map 
is  allowed  to  be  attenuated  by  natural  or  artificial,  consti¬ 
tuents  along  an  atmospheric  path.  In  addition,  radiation 
may  be  added  to  the  map  either  due  to  emission  or  scatter¬ 
ing.  The  natural  atmospheric  model  is  a  version  of  L0WTRAN6 
[1]  adapted  to  use  in  the  simulation.  The  standard  version 
of  L0WTRAN6  calculates  the  attenuation  caused  by  naturally 
occurring  molecular  and  aerosol  species.  It  also  calculates 
the  path  radiance  due  to  emission  from  each  species,  and 
aerosol  scattering  of  point  source  (eg.  solar)  radiation. 
It  is  assumed  that  natural  atmospheric  effects  will  be 
spatially  uniform  across  the  entire  image.  Further,  the 
spectral  variations  in  the  atmospheric  path  radiance  and 
attenuation  are  averaged  over  the  bandpass  of  interest  prior 
to  modifying  the  input  scene  radiance  map. 

The  battlefield  effects  model  is  a  modified  version  of  the 
ACTMAD  [2]  computer  code.  This  model  was  developed  by 
OptiMetrics  for  the  US  Army  Atmospheric  Sciences  Laboratory 
(ASL).  The  code  combines  the  ASL  smoke  model  ACT  II  [3] 
with  the  OptiMetrics  MADPUFF  [4,5]  smoke  model.  ACTMAD 
allows  up  to  6  varieties  of  smoke  munitions  and  12  total 
smoke  sources  to  be  considered  simultaneously.  The  sources 
may  be  positioned  in  any  configuration;  the  model  calculates 
a  spatial  map  of  the  obscurant  cloud  attenuation  and  path 
radiance.  Sky,  terrain,  and  point  radiation  sources  are 
considered  in  the  path  radiance  calculation.  A  notable 
feature  of  the  ACTMAD  model  is  that  it  includes  (as  a  user 
selected  option)  a  statistical  turbulent  smoke  transport  and 
diffusion  model.  Thus,  the  modeled  smoke  clouds  can  be  made 
to  appear  more  realistically  structured  than  is  possible 
with  the  more  conventional  gaussian-type  transport  and  dif¬ 
fusion  models. 

Because  we  have  calculated  the  optical  properties  of  the 
natural  and  battlefield  atmospheres  separately,  we  are  faced 
with  the  difficulty  of  combining  their  effects.  To  do  this. 


Figure  2-1.  Schematic  Representation  of  the  Simulation  Structure 


we  make  the  reasonable  assumptions  that  the  obscurant  cloud 
will  be  localized  along  the  atmospheric  path,  and  that  the 
natural  atmospheric  effect  within  an  obscurant  cloud  will  be 
negligible  compared  with  that  of  the  cloud  itself.  Based  on 
this,  the  atmospheric  path  is  divided  into  two  segments  with 
the  division  occurring  at  the  center  of  the  obscurant  cloud. 
The  L0WTRAN6  module  is  used  to  determine  the  optical  prop¬ 
erties  of  each  path.  The  net  path  optical  properties  are 
then  determined  according  to  the  following  equation: 


where 


L_(x,y)  -  L_  +  L  { x ,y )  t  +  l. 
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The  result  is  a  radiance  map  which  represents  the  input 
scene  as  it  would  appear  when  viewed  from  the  modeled  sensor 
location . 

The  sensor  model  basically  consists  of  two  parts.  A  simple 
optical  system  model  is  used  to  transform  the  apparent 
radiance  map  into  a  map  of  the  optical  power  (watts)  from 
each  scene  element  that  contributes  to  a  detector  signal. 
This  signal  map  is  then  Fourier  transformed  so  that  linear 
transfer  functions  representing  the  effect  of  sensor  subsys¬ 
tems  may  be  conveniently  applied.  White  noise  is  introduced 
following  the  sensor  detector,  and  is  then  filtered  by  sub¬ 
sequent  subsystems.  In  general,  the  sensor  model  utilizes 
the  same  subsystem  transfer  function  forms  as  the  NV&EOL 
Static  Performance  Model  for  Thermal  Viewing  Systems  [6]  . 
The  most  notable  exceptions  are  that  display  and  observer 


effects  are  simulated  rather  than  modeled  in  our  approach. 
Nevertheless,  the  sensor  descriptions  required  by  the  two 
models  are  nearly  identical  and,  for  simple  images  such  as 
Minimum  Resolvable  Temperature  (MRT)  bar  targets,  the  pre¬ 
dictions  of  the  two  models  should  be  comparable. 

A  special  capability  provided  in  the  model  is  the  ability  to 
perform  the  inverse  of  the  process  described  above.  In  this 
case,  a  measured  image  is  input  to  the  model  and  converted 
into  a  scene  radiance  map.  Basically  this  involves  applying 
the  inverse  of  the  measurement  sensor  transfer  functions, 
subtracting  the  atmospheric  path  radiance,  and  dividing  by 
the  atmospheric  path  transmission.*  Unfortunately,  the 
utility  of  this  method  is  limited  by  the  presence  of  noise 
in  the  measured  data,  the  accuracy  to  which  the  measurement 
sensor  can  be  represented,  and  the  ability  to  model  specific 
atmospheric  conditions. 


2 . 1  Atmospheric  Effects  Models 

Two  separate  atmospheric  effects  models  (L0WTRAN6  and 
ACTMAD)  are  incorporated  in  the  simulation  to  treat  natural 
and  battlefield  effects  respectively.  In  the  following  two 
subsections  we  will  provide  brief  descriptions  of  each 
model's  theoretical  basis;  however,  since  both  models  are 
fully  developed  entities  in  their  own  right,  the  reader  will 
be  referred  to  existing  technical  documentation  specific  to 
each  model.  Our  discussion  will  generally  be  limited  to 
details  peculiar  to  their  implementation  in  ISSM  as  well  as 
details  which  cannot  be  found  in  easily  obtainable  documen¬ 
tation  . 

2.1.1  Natural  ‘tmospheric  Effects  Model 

The  natural  atmospheric  effects  model  is  an  adaptation  of 
the  LOWTRAN 6  computer  code  [1]  which  was  developed  by  the  US 
Air  Force  Geophysics  Laboratory.  As  suggested  by  its  name, 
the  code  is  the  latest  version  in  a  series  of  atmospheric 
models  developed  over  the  past  several  years  [7,  8,  9,  10] . 

The  LOWTRAN  model  predicts  the  transmittance  and  path  radi¬ 
ance  for  an  arbitrary  path  geometry  at  a  spectral  resolution 
of  20  cra"l.  The  model  normally  represents  the  atmosphere 


*  Note:  In  the  current  ISSM  version  (1.0),  we  only  allow 
inversion  of  the  sensor  transfer  functions.  This  is  an 
oversight  which  will  be  corrected  in  version  2. 
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using  34  homogeneous  layers  which  describe  the  vertical  pro¬ 
files  of  temperature,  pressure,  water  vapor,  ozone,  and  at¬ 
mospheric  aerosols.  A  number  of  default  vertical  profiles 
are  stored  with  the  model;  the  user  may  alternatively  define 
a  specific  model  atmosphere  using  an  arbitrary  number  of 
layers.  In  addition  to  the  variable  atmospheric  quantities, 
LOWTRAN  also  treats  the  so-called  uniformly  mixed  gasses 
such  as  CO2  and  N20. 

Within  the  spectral  range  0.25-28.5  microns,  LOWTRAN  stores 
data  describing  the  optical  properties  of  all  the  naturally 
occur ing  gaseous  species,  and  a  large  number  of  aerosols 
including  rural  haze,  urban  haze,  fog,  rain,  snow,  tropo¬ 
spheric  aerosols,  and  volcanic  ash.  Using  this  optical 
data,  L0WTRAN6  calculates  the  attenuation  along  an  atmos¬ 
pheric  path  considering  the  following  factors: 

-  Molecular  absorption 

-  Molecular  scattering 

-  Aerosol  absorption 

-  Aerosol  scattering 

In  calculating  the  attenuation  due  to  scattering,  LOWTRAN 
assumes  that  the  signal  radiation  is  travelling  along  a 
single  ray  toward  a  point  receiver  with  an  infinitesimal 
f ield-of-view.  All  scattered  signal  radiation  is  considered 
"attenuated";  multiple  scattering  of  radiation  is  not  con¬ 
sidered.  Therefore,  the  theoretical  possibility  that  a 
scattering  media  may  have  a  non-uniform  Modulation  Transfer 
Function  ( MTF)  (that  is,  that  the  scattering  media  may  at¬ 
tenuate  signal  energy  at  different  spatial  frequencies  with 
different  efficiencies)  is  not  considered  by  LOWTRAN.  It 
should  be  noted,  however,  that  there  are  few  if  any  practi¬ 
cal  situations  involving  natural  atmospheric  effects  for 
which  the  above  consideration  is  of  any  significant  con¬ 
sequence  . 

LOWTRAN 6  calculates  the  atmospheric  path  radiance  con¬ 
sidering  the  following  factors: 

-  Molecular  emission  along  the  path  of  interest 

-  Aerosol  emission  along  the  path  of  interest 

-  Aerosol  single  scattering  of  plane-parallel  point 
sources  (eg.  solar  scattering) 


Specifically  excluded  from  consideration  are  multiply  scat¬ 
tered  photons,  atmospheric  emissions  scattered  into  the  path 
of  interest,  and  radiation  scattered  from  extended  sources 
such  as  the  ground.  Unfortunately,  at  the  longer  infrared 
wavelengths  of  interest,  the  later  two  effects  can  be  sig¬ 
nificant.  For  example,  in  the  8-12  micron  wavelength  region, 
the  emissions  from  the  sky  and  ground  hemispheres  can  easily 
dominate  the  contribution  from  the  sun.  Therefore,  L0WTRAN6 
may  significantly  under-estimate  the  atmospheric  path 
radiance.  As  is  discussed  in  section  2.2,  this  deficiency 
is  of  no  consequence  in  the  current  version  of  the  ISSM 
since  only  spatially  uniform  atmospheric  emissions  are  con¬ 
sidered;  the  linear  sensor  modeling  approach  does  not  con¬ 
sider  signal  processing  dynamic  range  limitations  so  that 
the  degradation  caused  by  spatially  uniform  emissions  can 
always  be  eliminated  by  the  (modeled)  sensor.  However,  if 
later  versions  of  the  ISSM  consider  non-uniform  atmospheric 
effects  (eg.  consideration  of  the  horizon),  then  the  defic¬ 
iency  should  be  corrected  using  the  ACTMAD  radiative  trans¬ 
fer  modeling  techniques  (which  include  the  effect  of  ex¬ 
tended  sources)  discussed  in  the  next  section. 

2.1.2  Battlefield  Effects  Model 

The  ACTMAD  [2]  battlefield  effects  model  is  a  derivative  of 
the  ACT  II  [3]  smoke  obscuration  model  developed  by  the  US 
Army  Atmospheric  Sciences  Laboratory  (ASL) .  The  ACT  II  model 
treats  the  transport  and  diffusion  of  obscurant  clouds  in  a 
deterministic  manner  using  a  gaussian  plume  methodology.  In 
order  to  model  the  spatial  and  temporal  inhomogeneities 
present  in  obscurant  clouds  generated  under  realistic  field 
conditions,  OptiMetrics  developed  the  Microscale  Atmospheric 
Diffusion  of  Puffs,  or  MADPUFF  methodology.  This  methodology 
has  been  incorporated  in  the  ACT  II  computer  code  as  a  user 
selected  option;  the  resulting  ACTMAD  computer  code  was 
developed  by  OptiMetrics  for  ASL. 

2. 1.2.1  ACT  II  Smoke  Obscuration  Model 

The  ACT  II  smoke  obscuration  effects  model  [3]  allows  for 
the  calculation  of  transmission  and  path  radiance  along  an 
arbitrary  path  through  an  obscurant  cloud.  Obscurant  sources 
are  represented  using  a  time  series  of  puffs  (usually  pro¬ 
duced  at  a  rate  of  one  per  second)  with  a  spatially  gaussian 
mass  distribution.  The  source  characteristics  are  deter¬ 
mined  by  the  initial  size,  mass,  temperature,  and  composi¬ 
tion  of  each  puff.  The  transport  and  diffusion  of  each  puff 
is  modeled  using  an  empirical  gaussian  plume  methodology 
developed  by  Pasquill  [11].  In  this  method,  each  puff  is 
allowed  to  travel  with  the  mean  atmospheric  flow  and  diffuse 
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according  to  empirical  laws  driven  by  the  atmospheric  sta¬ 
bility.  Puffs  are  also  allowed  to  rise  based  on  the  buoyant 
force  afforded  by  their  temperature.  The  ACT  II  model  treats 
the  transport  of  all  puffs  identically  so  that  the  resulting 
obscurant  cloud  structure  is  very  smooth  and  regular;  the 
only  cloud  structure  arises  from  modeled  temporal  variations 
in  the  obscurant  source. 

Once  a  particular  obscurant  cloud  geometry  is  realized  using 
the  transport  and  diffusion  models,  the  transmission  and 
path  radiance  may  be  calculated  using  a  very  complete  radia¬ 
tive  transfer  model.  The  optical  properties  of  the  obscur¬ 
ant  (extinction  coefficient,  single  scattering  albedo,  and 
scattering  phase  function)  must  be  input  by  the  user.  Calcu¬ 
lations  are  performed  using  only  one  set  of  optical  para¬ 

meters;  the  selection  of  these  parameters  determines  whether 
the  calculations  are  valid  at  a  single  wavelength,  or 
whether  they  represent  a  broadband  average  result.  The 

scope  of  the  radiative  transfer  model  is  illustrated  in 
Figure  2-2.  Transmission  calculations  are  performed  account¬ 
ing  for  aerosol  absorption  and  scattering;  as  is  the  case 
with  the  LOWTRAN  calculations,  effects  leading  to  a  non- 
uniform  MTF  such  as  multiple  scattering  are  not  accounted 

for.  The  path  radiance  calculations  are  considerably  more 
extensive  than  those  of  LOWTRAN.  ACT  II  considers  the 

complete  single  scattering  problem  and  accounts  for: 

-  Radiation  emitted  along  the  path  of  interest; 

Smoke  emissions  scattered  into  the  path  of 
in terest ; 

Plane  parallel  source  emission  scattered  into  the 
path  (eg.  solar  scattering);  and 

Radiation  from  the  sky-ground  sphere  scattered  into 
the  path. 

To  account  for  variations  in  the  sky  and  ground  radiation,  a 
sectoring  scheme  is  used.  The  user  may  divide  the  sky  into 
any  number  of  sectors,  and  supply  the  diffuse  radiance  con¬ 
tribution  from  each  sector;  the  ground  is  represented  by  a 
homogeneous  grey  body. 

An  addition  to  the  ACT  II  radiative  transfer  model  is  a  very 
simple  path  radiance  approximation.  The  purpose  of  this 
approximation  is  to  allow  very  rapid  calculations  while 
still  including  terms  indicative  of  the  actual  obscurant 
emission  and  scattering  terms.  The  radiance  equation  is: 


Figure  2-2.  Illustration  Showing  the  Effects  Considered  in  the 
Radiative  Transfer  Model 


(2) 


Lpu  y)  =  (i*°“Ts(x»y)/  [dgP  •>  ♦  lbb(ts)  (1— )J 


where  u>  =  Obscurant  single  scattering  albedo 

I  ■  Total  sky/ground  irradiance  at  the  obscur- 
3  ant  cloud  location  (W/m2) 

P  *  Scattering  phase  function  weighted  by  the 

radiance  distribution  in  the  sky  and  ground 
hemispheres  (Sr-1) 

T  »  Average  obscurant  temperature  input  by  the 
3  user 

Lno  »  Black  body  radiance  function  (W/m2-Sr) 


The  approximation  is  most  useful  when  initially  setting  up  a 
problem  for  consideration,  and  representative  results  are 
desired  quickly.  The  full  single  scattering  model  should  be 
used  when  actual  simulations  are  attempted. 

2. 1.2. 2.  MADPUFF  Transport  and  Diffusion  Methodology 

The  MADPUFF  method  is  described  in  several  papers  and  Opti- 
Metrics  reports  [2,4,5].  The  method  is  based  on  work  per¬ 
formed  by  S.R.  Hanna  [12,13].  The  modeling  approach  assumes 
that  the  velocity  (u)  of  a  fluid  parcel  at  point  (x,y,z)  at 
time  t  can  be  represented  as  the  sum  of  a  mean  wind  compo¬ 
nent  U  and  a  turbulent  fluctuation  component  u': 


u(x,y,z,t)  *  U(x,y, z, t)  +  u’(x.y,z,t) 


The  turbulent  fluctuation  component  u'  is  also  assumed  to  be 
the  sum  of  two  components  such  that: 


u'  (t  +  t)  »  u'(t)  R(  t)  +  u* 


where  R( x)  ■  Lagrangian  autocorrelation  function  for 

time  lag  x 

■  exp  [-x/Tl],  Tr«  Lagrangian  time  scale 


u"  *  normally  distributed  randan  turbulence 

component  with  standard  deviation  given  by: 

v-  v 11 "  r2( 

Hanna  combines  theory  and  empirical  findings  [14]  to  relate 
the  Lagrangian  time  scale  to  more  readily  observed  meteo¬ 
rological  quantities: 

u:  Z  <  Zt 

(5) 

v:  Z  <  Z^ 

w:  100  m  <  Z  <  Z^ 

T.  -  0.42  (Z/o  )  (u/u*)  for  (w:  Z  <  100  m) 

Lw  w 

where 

Z^  *  mixing  depth 

u  »  mean  wind 

°u  v  w  *  Lagrangian  standard  deviation  for  u,v,w 

'  '  components,  respectively 

u*  »  friction  velocity 

Z  *  height  above  ground 


In  the  MADPUFF  method,  gaussian  smoke  puffs  are  advected  by 
the  randomly  fluctuating  wind  vector  u;  the  fluctuating  part 
of  the  wind  (u”)  is  simulated  using  a  time  series  of  random 
numbers.  Note  that  although  the  puffs  are  advected  randomly, 
their  diffusion  is  still  governed  by  the  same  empirical  re¬ 
lations  used  in  ACT  II.  However,  since  the  puff  diffusion 
model  used  in  ACT  II  is  intended  to  represent  both  the  cloud 
meandering  as  well  as  the  puff  growth,  the  puff  diffusion 
terms  used  with  MADPUFF  are  smaller  by  a  factor  of  1/3  [11] 
since  the  random  puff  motion  is  considered  separately. 

The  only  remaining  difference  between  the  ACT  II  mode  and 
the  MADPUFF  mode  is  the  puff  buoyancy  model.  OptiMetrics 
developed  an  empirical  buoyancy  model  based  on  data  obtained 
during  the  Smoke  Week  IV  and  Jefferson  Proving  Ground  trials 
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[2];  this  model  is  used  when  the  MADPUFF  mode  is  selected. 
The  original  ACT  II  buoyancy  model  is  used  when  the  ACT  II 
mode  is  selected. 

2.2  Sensor  Model 

Since  the  sensor  model  was  developed  for  this  simulation, 
the  discussion  will  of  necessity  be  more  detailed  than  the 
previous  sections.  The  ISSM  sensor  model  structure  was  de¬ 
picted  in  Figure  2-1.  A  simple  optical  system  model  is  used 
to  determine  the  efficiency  with  which  scene  radiation  is 
directed  to  the  system  detector(s).  The  resulting  signal 
map  is  Fourier  transformed  so  that  linear  transfer  function 
representations  of  the  sensor  subsystems  may  be  easily  ap¬ 
plied.  For  convenience,  most  of  the  transfer  functions 
incorporated  in  the  model  are  identical  to  those  used  in  the 
NV&EOL  Static  Performance  Model  [6] .  This  means  that  sensor 
descriptions  used  with  the  NVL&EOL  model  will  also  be  ap¬ 
plicable  for  use  in  the  simulation.  White  noise  (although 
any  character  noise  may  be  accommodated  by  the  model)  is 
inserted  following  the  system  detector;  the  detector  is  as¬ 
sumed  to  be  the  sensor's  limiting  noise  source.  After  all 
the  transfer  functions  have  been  applied,  a  sensor  output 
image  is  formed  using  an  inverse  Fourier  transform.  This 
image  may  be  displayed  and  evaluated  by  the  user. 

2.2.1  Optical  System  Model 

The  optical  system  is  characterized  in  terms  of  a  simple 
lens  defined  by  its  f-number,  focal  length,  and  transmis¬ 
sion,  and  a  filter  defined  by  its  spectral  transmission 
within  the  sensor  bandpass.  When  combined  with  the 

characteristics  of  a  single  detector  element,  the  signal 
from  each  scene  pixel  in  units  of  watts  is  given  by: 

"A.  /  X 

S(x,y)  -  - 2 - 5-  /  *  R(X)  xn (  A)  L.(x,y)  dX 

4  (f/#T  J  X  U  A  (6) 


S(x,y) 

* 

Detector 

Signal  Map  (watts) 

Ad 

at 

Detector 

Area  (m^) 

f/# 

- 

Optical 

System  F-number 

R<  X) 

- 

De tec tor 

Normalized  Spectral  Response 

. 

Optical 

System  Spectral  Transmission 

L»(x,y)  *  Apparent  Scene  Radiance  Map 

A  <w/m2-sr-micron) 


and  in  the  current  version  of  ISSM,  the  scene  radiance  map 
is  assumed  to  be  spectrally  uniform. 


2.2.2  2-Dimensional  Image  Fourier  Transform 

The  standard  Numerical  Analysis  Library  (NAL)  [15]  routines 
SFFT2A  and  SFFTA  are  used  to  perform  the  Fast  Fourier  Trans¬ 
form  (FFT).  Both  routines  perform  the  FFT  on  n-dimensional 
single  precision  complex  arrays;  the  difference  between  the 
two  routines  is  that  SFFT2A  requires  the  number  of  elements 
in  each  array  dimension  to  be  a  power  of  two  while  no  such 
restriction  exists  for  SFFTA. 

Before  the  signal  array  S  can  be  transformed,  the  data  must 
be  made  periodic  so  that  spurious  "edge  effects"  are  not 
introduced.  This  is  done  by  reflecting  the  data  as  illus¬ 
trated  in  Figure  2-3.  Notice  that  this  procedure  does  not 
introduce  spurious  information  into  the  data;  because  of  the 
symmetry,  all  odd  Fourier  components  in  both  dimensions  will 
be  zero.  Of  course,  one  may  take  advantage  of  this  by  per¬ 
forming  computations  using  only  the  non-zero  (even)  Fourier 
components. 

When  performing  the  Fourier  Transform  of  a  2-dimensional 
complex  array  of  size  nxm,  the  result  is  a  complex  array  of 
the  same  size.  However,  when  the  original  array  is  real, 
then  the  transformed  array  will  have  conjugate  symmetry 
about  each  dimension.  Thus,  all  the  data  information  is 
really  contained  in  a  (n/2)x(m/2)  complex  array;  computa¬ 
tions  need  only  be  performed  using  one  quarter  of  the 
transformed  array.  This  means  that  reflecting  the  image 
data  (Figure  2-3)  has  no  effect  on  the  number  of  compu¬ 
tations  which  must  be  performed  in  the  Fourier  domain. 

Finally,  it  should  be  noted  that  techniques  have  been 
developed  to  allow  real  arrays  of  size  n  to  be  transformed 
using  a  complex  array  of  size  n/2  [16].  Of  course,  this 

technique  can  eliminate  the  computational  penalty  involved 
in  taking  the  Fourier  Transform  of  the  reflected  image.  We 
have  not  implemented  this  approach  in  the  ISSM;  however,  if 
very  large  images  (eg.  512x512)  are  ultimately  processed, 
then  consideration  of  this  technique  may  prove  to  be 
important. 


(1,1) 

(n,  1) 

(n,  1) 

(1,1) 

ORIGINAL 

IMAGE 

(1,  m) 

(n,m) 

(n,m) 

(l,m) 

(l,m) 

(n,m) 

(n,m) 

(l,m) 

(1,1) 

(n,  1) 

(n,  1) 

(1,1) 

Figure  2-3.  Illustration  Showing  the  Image  Reflection 
Scheme  used  to  Eliminate  Edge  Effects  in 
the  Fourier  Transform  Operation 


2.2.3.  Sensor  Subsystem  Transfer  Functions 

Each  sensor  subsystem  (optics,  detector,  etc.)  is  repre¬ 
sented  by  a  real  or  complex  transfer  function.  Because  we 
are  operating  in  the  frequency  domain,  the  linear  effects  of 
each  system  element  are  represented  by  the  product  of  the 
transformed  image  and  a  transfer  function.  The  complete  set 
of  pre-defined  subsystem  transfer  functions  is  listed  in 
Table  2-1.  Not  mentioned  in  the  table  is  the  fact  that  for 
several  subsystems,  the  user  has  the  option  of  numerically 
defining  the  transfer  function.  We  will  briefly  comment  on 
each  transfer  function  in  the  following  subsections. 

2. 2. 3.1.  Optical  System  Transfer  Functions 

Two  optical  system  transfer  function  forms  are  provided 
which  describe  diffraction  by  a  circular  aperture  and  blur 
due  to  optical  system  aberrations  [6] .  A  gaussian  transfer 
function  form  is  used  to  represent  the  blur;  it  is  specified 
by  the  1/e  half-width  of  the  corresponding  gaussian  point 
spread  function.  Note  that  in  many  real  cases,  this  is  not 
a  very  realistic  representation  of  optical  system  aberra¬ 
tions.  For  these  cases,  the  user  may  provide  a  numerical 
optical  system  transfer  function;  transfer  functions  cor¬ 
responding  to  many  common  aberrations  are  found  in  refer¬ 
ences  17  and  18. 

2. 2. 3. 2.  Atmospheric  Turbulence  Transfer  Function 

The  atmospheric  turbulence  transfer  function  used  in  the 
model  originated  with  Fried  [19,20],  We  have  chosen  the 
near-field  formulation  of  the  short  exposure  case.  The 
near-field  assumption  encompasses  most  army  scenarios  of 
interest,  and  the  short  exposure  assumption  is  consistent 
with  the  rapid  frame  rates  of  most  imaging  systems. 

Representing  the  effect  of  atmospheric  turbulence  using  a 
Modulation  Transfer  Function  (MTF )  is  at  best  a  rather  crude 
approximation  for  several  reasons.  First,  turbulence  is 
dynamic  so  that  it  will  affect  a  particular  image  feature 
differently  at  each  instant  of  time.  For  example,  there  is 
a  finite  probability  of  obtaining  a  diffraction-limited  view 
of  a  particular  image  feature;  the  MTF  we  use  is  really  an 
average  of  the  degradation  caused  in  an  ensemble  of  short 
exposure  images.  The  second  point  is  that  turbulence  effects 
are  not  spatially  symmetric  so  that  a  very  specific  Complex 
or  Optical  Transfer  Function  (OTF)  rather  than  an  MTF  is 
required  even  at  a  single  instant  of  time.  Indeed,  a  single 
MTF  description  is  strictly  valid  only  over  a  small  portion 
of  each  image  called  the  isoplanatic  patch  [21]  .  However, 
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Table  2-1.  Sensor  Subsystem  Transfer  Function  Definitions 
(continued) 


Atmospheric  Turbu lence 


,  Dkxy  \5/3 
H(Rx,ky)  -  e  3-44[™*\/ 


0_0  YXPkxy\ 

.0  o*51  1000  J 


Rq  *  atmospheric  coherence  length  (m) 
-  fo.159  C2/,  2J  .■  V  rl  ~0'6 


L  "V1*"  V  J 

D  *  optical  system  aperture  diameter  (m) 

C2  =  refractive  index  structure  parameter  (m“^/^) 
n 

r  *  sensor-target  range  (m) 

Detector  -  Temporal  Transfer  Function 


H(kx)  » 


1.0  +  j 


kxvN 


FOV  FOV  F  n 
x  y  r  ov 


n  I FOV  n 
p  y  sc 


v  =  scan  velocity  (mrad/s) 

?x  =  sensor  along-scan  f ield-of-view  (mrad) 
/y  =  sensor  cross-scan  f ield-of-view  (mrad) 

*  frame  rate  (s-*) 

*  overscan  ratio 

/ 

*  number  of  detectors  in  parallel 


IFOVy  ■  sensor  instantaneous  f ield-of-view  in  the 
y  cross-scan  direction  (mrad) 


scan  efficiency 


Table  2-1.  Sensor  Subsystem  Transfer  Function  Definitions 
(continued) 


Detector-Spatial  Transfer  Punction 


H  lkx  'ky  )  ' 

sinc{*  IPOV  k  J  sinc(*  IFOV  k  ) 
xx  y  y 

IFOVx 

sensor  instantaneous  f ield-of-view  in  the 
along-scan  direction  (mrad) 

sinc(x)  * 

sin(x)/x 

Electronics 

h(kx) 


uoT77£I\  i  +  j 


(t) 


electronics  cut-off  frequency  (Hz) 
electronics  cut-on  frequency  (Hz) 


Electronic  Boost 


H(kxJ  -  1.0  +  (B^6°-~  1*0  -  cos  ( 

\  max  i 


frequency  of  maximum  boost  (Hz) 


boost  at  frequency  f 


Table  2-1.  Sensor  Subsystem  Transfer  Function  Definitions 
(continued) 

LEDs 

H  (k  ,k  )  *  sine  (if  LED  k  )  sinc(*  LED  k  ) 
x  y  xx  y  y 

LEDX  =  along-scan  LED  dimension  (mrad) 

LEDy  *  cross-scan  LED  dimension  (mrad) 


f  *  Nyquist  frequency  determined  by  the  sampling 
nyq  rate  (Hz) 

Stabilization 

.ivy  •  W. -(v>‘) 

Sx  *  along-scan  stabilization  parameter  (mrad2) 

Sy  »  cross-scan  stabilization  parameter  (mrad2) 


if  we  were  to  accumulate  the  specific  OTFs  which  accurately 
describe  the  effect  of  turbulence  on  short  exposure  images, 
we  would  find  the  all  portions  of  the  image  are,  on  the 
average,  affected  in  a  way  described  by  the  MTF  listed  in 
Table  2-1. 

Finally,  a  few  of  the  parameters  on  which  the  turbulence 
transfer  function  depends  may  be  unfamiliar  to  the  reader. 
The  coherence  length  is  the  lateral  distance  over  which  a 
wavefront  remains  coherent  after  passing  through  the  turbu¬ 
lent  media.  An  interesting  interpretation  of  the  coherence 
length  is  that  the  attainable  resolution  in  a  long  exposure 
image  is  approximately  equal  to  that  of  a  diffraction  limit¬ 
ed  system  whose  aperture  diameter  is  equal  to  the  coherence 
length.  A  formula  for  calculating  the  coherence  length  af¬ 
forded  by  a  homogeneous  atmospheric  path  is  listed  in  Table 
2-1.  The  descriptor  of  the  turbulence  strength  is  the 
refractive  index  structure  parameter  (c  2).  Reference  22 
provides  an  excellent  discussion  of  turbulence  in  the 
atmosphere  including  the  typical  diurnal  behavior  of  C  1  . 
For  convenience,  some  typical  C  2  values  are  listed  bel3w: 


Time  of  Day 


Cn2  (m~2/3) 


sunrise  or  sunset 

midday 

midnight 

extreme  midday  condition 


1. E-15  or  lower 

2. E-14  -  6.E-14 

1. E-14  -  2.E-14 

2. E-13  -  5.E-13 


2. 2. 3. 3.  Detector  Transfer  Functions 


Two  transfer  functions  have  been  included  to  represent  the 
temporal  and  spatial  detector  response.  The  temporal  re¬ 
sponse  is  modeled  by  a  complex  transfer  function  represent¬ 
ing  a  single  pole  low-pass  filter.  The  spatial  response  of 
the  detector  is  modeled  using  the  2-dimensional  MTF  of  a 
rectangular  aperture.  Of  course,  this  means  that  the  de¬ 
tector  elements  are  assumed  to  respond  uniformly  over  their 
entire  area. 


It  should  be  noted  that  in  this  version  of  the  ISSM,  we  have 
not  considered  the  effect  of  cross-scan  sampling.  In  most 
thermal  imaging  sensors,  cross-scan  sampling  occurs  because 
the  image  is  scanned  in  a  raster  pattern  by  one  or  more  dis¬ 
crete  detectors.  Although  the  effect  will  be  considered 
along  with  other  non-linear  phenomena  in  later  versions  of 


ISSM,  it  would  be  possible  to  implement  an  analogous  pro¬ 
cedure  to  that  used  to  represent  along-scan  sampling  (sec¬ 
tion  2. 2. 3. 7)  should  this  be  desired  by  the  reader. 

2. 2. 3. 4.  Electronics  Transfer  Functions 

Two  transfer  function  forms  are  included  to  represent  the 
sensor  electronics;  a  complex  transfer  function  representing 
single  pole  low-pass  and  high-pass  filters,  and  a  MTF  repre¬ 
senting  a  high  frequency  boost  function.  The  form  of  the 
boost  function  is  exactly  that  used  by  NV&EOL  [6] . 

2. 2. 3. 5.  Light  Emitting  Diode  (LED)  Transfer  Functions 

Virtually  all  existing  thermal  imaging  systems  which  incor¬ 
porate  parallel  detector  arrays  use  LEDs  to  either  display 
the  image  directly,  or  as  part  of  an  Electro-Optic  (EO) 
multiplexer.  The  LEDs  are  assumed  to  be  rectangular;  their 
spatial  transfer  function  is  identical  in  form  to  that  of 
the  sensor  detectors.  No  temporal  filtering  by  the  LEDs  in 
the  along-scan  dimension  is  considered,  and  any  sampling 
effects  which  involve  the  LEDs  are  ignored. 

2. 2. 3. 6.  Vidicon  Transfer  Function 

No  default  transfer  function  form  is  included  for  a  vidicon; 
the  user  must  input  a  numerical  definition.  The  vidicon 
(which  is  the  video  multiplexer  in  most  current  technology 
infrared  imaging  systems)  is  really  a  complete  EO  system  in 
its  own  right.  Typically,  it  will  have  an  asymmetric  trans¬ 
fer  function  that  may  include  boost  in  the  along-scan  dimen¬ 
sion.  Thus,  it  is  difficult  to  provide  a  suitable  default 
transfer  function  form. 

2. 2. 3. 7.  Digital  Multiplexer  (Sampling)  Transfer  Function 

In  a  digital  multiplexer  system,  the  LEDs  and  vidicon  which 
make  up  the  EO  multiplexer  are  replaced  by  an  electronic 
circuit  which  combines  the  output  of  each  detector  in  a 
single  video  stream.  In  a  typical  electronic  multiplexer, 
each  detector  signal  is  applied  to  a  separate  charge  accum¬ 
ulator,  and  its  signal  is  integrated  for  one  sample  period. 
The  accumulated  charges  are  then  gated  to  a  serial  charge 
transport  register.  By  reading  the  entire  contents  of  this 
register  within  one  sample  period,  a  serial  video  stream  is 
formed.  Note  that  the  video  is  serial  in  the  cross-scan 
dimension;  thus  the  process  introduces  sampling  in  the 
along-scan  direction. 


The  transfer  function  given  in  Table  2-1  accounts  for  both 
the  integration  and  sampling  of  the  assumed  electronic 
multiplexer.  The  sine  function  arises  from  assuming  that 
the  charge  accumulator  is  represented  by  an  ideal  or  "box¬ 
car"  integrator.  Finally,  the  Nyquist  frequency  is  simply 
one  half  the  sampling  rate;  the  transfer  function  is  written 
assuming  that  the  sample  period  and  integration  period  are 
exactly  matched. 

2. 2. 3. 8.  Stabilization  Transfer  Function 


The  stabilization  transfer  function  may  be  defined  either 
numerically,  or  using  the  default  gaussian  form.  The  input 
parameters  for  the  gaussian  form  are  exactly  the  same  as 
those  used  in  the  NV&EOL  model  [ 6 ]  .  Note  that  the  para¬ 
meters  Sx  and  Sv  may  be  related  to  the  stabilization-limited 
point  spread  function  1/e  halfwidth  as  indicated  in  the 
optical  system  blur  transfer  function  formula. 


2.2.4  Detector  Noise  Model  and  Implementation 


The  ISSM  allows  noise  to  be  added  at  a  point  just  following 
the  system  detector(s).  The  technique  used  is  to  compute  a 
random  sample  of  noise  and  then  add  this  sample  to  the  scene 
signal  in  the  frequency  domain.  All  sensor  subsystems  fol¬ 
lowing  this  point  of  noise  insertion  then  act  on  both  the 
signal  and  noise.  Adding  the  noise  in  the  frequency  domain 
is  very  convenient  in  that  it  is  trivial  to  account  for  any 
arbitrary  noise  power  spectrum,  and  additional  Fourier 
Transform  operations  are  avoided.  However,  to  avoid  confu¬ 
sion,  it  should  be  noted  that  the  current  ISSM  configuration 
uses  a  simple  white  noise  model. 


Use  of  a  white  noise  model  is  justified  in  many  cases  of 
interest.  For  example,  photovoltaic  or  photoconduct ive 
detectors  are  well  represented  by  a  white  noise  spectrum 
over  a  broad  range  of  intermediate  frequencies  limited  only 
by  the  detector  temporal  response  at  very  high  frequencies, 
and  1/f  noise  at  very  low  frequencies.  In  practical  imaging 
sensor  applications,  the  system  electronics  usually  restrict 
the  system  bandwidth  well  below  the  detector  high  frequency 
response.  Further,  the  1/f  noise  is  usually  filtered  by  ac 
coupling  between  the  detectors  and  their  preamplifiers. 

The  modeled  noise  power  spectral  density  is  given  as 
follows: 


A 


(7) 


where 


*n  ■  noise  power  spectral  density  (W2/Hz) 

A^  *  detector  element  area  (cm2) 

n  ■  number  of  detectors  in  series 
s 

F  *  frame  rate  (s“*) 

r 

*  eye  integration  time  (s) 

D*  *  detector  peak  detectivity  (cm-( Hz ) 1//2/w) 


There  are  several  points  concerning  this  expression  which 
deserve  comment.  First,  it  should  be  emphasized  that  the 
detector  D*  characteristic  is  spectrally  dependent.  It  is 
important  that  the  noise  be  calculated  using  the  D*  relevant 
at  the  point  of  peak  spectral  response  (as  defined  by  the 
response  characteristic  R(  X)  in  equation  6)  in  order  to  be 
consistent  with  the  way  that  the  signal  is  defined.  Second, 
the  quantity  ns  is  present  to  account  for  Time  Delay  and 
Integration  (TDI )  techniques  which  can  be  used  to  reduce 
noise  in  sensors  which  have  multiple  detectors  positioned  in 
series.  Finally,  the  term  (FrTe)  attempts  to  correct  for 
the  fact  that  the  simulation  can  only  compute  a  single, 
static  image  for  display  to  an  observer.  In  a  real  sensor, 
a  number  of  image  frames  with  statistically  independent 
noise  samples  may  be  presented  to  the  observer  within  one 
eye  integration  time.  Because  of  this,  the  observer  using 
the  real  sensor  can  obtain  a  better  image  signal-to-noise 
ratio  than  can  be  obtained  with  the  static  image.  Of  course, 
the  approximate  correction  for  this  is  the  term  given  above. 

The  actual  sample  of  noise  which  is  added  to  an  image  pixel 
in  the  Fourier  doman  is  computed  as  follows: 


Noise^ 


(V)  V  *  (V) 


Nk+1 


where  Noise,  *  noise  added  at  a  point  in  the  fourier 

domain  (W/fHj)1'2) 

Nj^  *  normally  distributed  random  number 
with  zero  mean  and  unit  variance 
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We  are  simply  computing  a  field  of  normally  distributed 
complex  random  numbers  whose  mean  is  zero  and  whose  variance 
is  the  white  noise  power  spectral  density  ,  Note  that  it 
would  be  a  trivial  matter  to  include  a  non-white  power 
spectrum  in  equation  8  should  that  become  necessary  in  the 
future. 


3.0.  User's  Guide 


The  routines  which  comprise  the  Imaging  Sensor  Simulation 
Model  (ISSM)  may  be  grouped  into  six  general  categories: 

1.  control  routines, 

2.  frame  I/O  interface  routines, 

3.  sensor  module  routines, 

4.  atmospheric  effects  module  routines, 

5.  battlefield  effects  module  routines,  and 

6.  image  display  routines. 

Tables  3-1  through  3-6  provide  descriptions  of  each  subpro¬ 
gram  included  in  each  of  these  six  categories. 

3.1.  Using  the  Simulation 

3.1.1.  Logical  Unit  Assignments 

A  large  number  of  auxiliary  input/output  unit  assignments 
are  necessary  to  use  all  the  features  of  the  simulation. 
For  example,  data  libraries  may  be  used  to  input  complicated 
sets  of  input  parameters  such  as  those  required  by  the  sen¬ 
sor  module  and  the  battlefield  effects  module.  Similarly, 
the  image  output  may  take  a  number  of  forms  such  as  video 
display,  printer  plots,  and  formatted  or  unformatted  files 
on  tape . 

Table  3-7  presents  a  list  of  the  input/output  units  cur¬ 
rently  used  in  the  simulation.  A  brief  description  of  the 
unit's  purpose  and  the  form  of  the  input  or  output  is  pro¬ 
vided. 

3.1.2.  Menu  Hierarchy 

The  simulation  is  structured  around  a  menu  tree  which  allows 
the  user  to  access  the  various  modules  in  any  logical  order. 
Each  menu  returns  control  to  the  menu  from  which  it  was 
called  after  it  either  performs  its  functions  or  is  explic¬ 
itly  told  to  return  by  the  user.  Since  menus  may  be  nested 
several  levels  deep,  it  is  important  for  the  user  to  become 
familiar  with  the  menu  hierarchy.  Figures  3-1  through  3-5 
provide  a  schematic  representation  of  the  menu  hierarchy  in 
the  simulation  which  will  enable  the  user  to  visualize  how 
one  can  pass  through  the  menus  to  access  certain  features  of 
the  model.  Figure  3-1  provides  the  main  menu  structure  and 
the  menus  immediately  accessed  from  it.  Figures  3-2  through 
3-5  detail  the  menu  structure  of  the  four  major  modules: 
the  frame  I/O  interface,  sensor,  atmospheric  effects,  and 
battlefield  effects  modules. 


Table  3-1.  Control  Routines  Included  in  Simulation 


These  routines  set  up  the  initial  conditions  and  provide  an 
interface  between  the  main  menu  and  individual  modules. 


Routine 


Description 


MAIN 


Reads  in  sensor  library,  sets  terminal 
type,  selects  main  menu  options. 


BLOCK  DATA 


Default  values  for  variables  input  by 
menus. 


LIBR 


Interface  with  frame  I/O  routines, 


SNSR 


Interface  with  sensor  routines. 


ATMOS 


Interface  with  atmospheric  effects 
module. 


BATTLE 


Interface  with  battlefield  effects 
module. 


PRLINE 


Interface  with  image  display  routines 


Routine 


Description 


LIBR 


FRREAD 


Selects  options  for  input  or  output  to  a 
data  library. 


Reads  frame  from  file  and/or  selects  a 
part  of  frame  for  processing. 


TAPEO 


Outputs  frame  to  a  library 


Table  3-3.  Routines  Included  in  Sensor  Effects  Module 


Routine 


SNSR 

SENCHR 

TFSWCH 

TRNSFM 

I 

OTFS 

HOD  IF 

HBLUR 

HTURB 

HDET1 

ZDET2 

ZELEC 

HBOOST 

ZSAMP 

SRESP 

SFFTA 

SFFT2A 

ALINEY 


Description 

Displays  sensor  menu,  calls  routines  for 
appropriate  options. 

Modifies  sensor  characteristics  according 
to  user  inputs. 

Selects  which  transfer  functions  will  be 
applied  to  represent  sensor  effects. 

Calls  Fourier  transform  routines,  optical 
transfer  function  routine. 

Sets  up  and  calls  optical  transfer 
functions. 

Diffraction  optical  transfer  function. 

Blur  optical  transfer  function. 

Atmospheric  turbulence  transfer  function. 

Spatial  detector  transfer  function. 

Temporal  detector  transfer  function. 

Electronics  transfer  function. 

Electronic  boost  transfer  function. 

Along-scan  sampling  transfer  function. 

Calculates  optical  system  response. 

Fast  Fourier  transform  for  an  arbitrarily 
dimensioned  array. 

Fast  Fourier  transform  for  an  array  whose 
dimensions  are  powers  of  2. 

Linear  interpolation  routine. 

Integration  routine. 


QTFE 


Routine 


Description 


LWTRN6  Calls  subroutines  to  compute  transmit¬ 

tance  and  radiance. 

NSMDL  For  user-defined  atmospheric  or  aerosol 

profiles. 

SUBSOL  Calculates  the  subsolar  point  angles 

based  upon  time  and  day. 

STDMDL  Sets  up  atmospheric  profiles  of  attenu¬ 

ator  densities. 

AERPRF  Computes  scaling  factor  profiles  for 

aerosols. 

GEO  Driver  for  air  mass  subroutines.  Calcu 

lates  attenuator  amounts  for  the  slant 
path. 

SSGEO  Obtains  attenuator  amounts  from  scatter 

ing  points  along  optical  path  to  the 
extraterrestrial  source. 

EXABIN  Loads  aerosol  extinction  and  absorption 

coefficients  for  the  appropriate  model 
and  relative  humidity. 

TRANS  Calculates  transmittance,  atmospheric 

radiance,  and  solar/lunar  scattered 
radiance  for  slant  path. 


(continued) 


Routine 


Description 


GEO 

GEOINP 

FNDHMN 

REDUCE 

FDBETA 

RFPATH 

FILL 

LAYER 

RADREF 

FINDSH 

i 

SCALHT 

t 

|  ANDEX 


Driver  for  air  mass  subroutines.  Calcu¬ 
lates  attenuator  amounts  for  the  slant 
path. 

Interprets  geometry  input  parameters  into 
the  standard  form  HI,  H2,  ANGLE,  and  LEN. 

Calculates  HMIN,  the  minimum  altitude 
along  the  path  and  PHI,  the  zenith  angle 
at  H2. 

Eliminates  slant  path  segments  that  ex¬ 
tend  beyond  the  highest  profile  altitude. 

Calculates  angle,  given  HI,  H2,  and  BETA 
by  iteration. 

Determines  the  refracted  path  and  the 
absorber  amounts  through  all  the  layers. 

Defines  the  boundaries  of  the  slant  path 
and  interpolates  densities  at  these 
boundaries . 

Calculates  the  path  and  amounts  through 
one  layer. 

Computes  radius  of  curvature  of  the 
refracted  ray  for  a  horizontal  path. 

Finds  layer  boundaries  and  scale  height 
at  ground  for  index  of  refraction. 

Calculates  scale  height  of  index  of 
refraction. 

Computes  index  of  refraction  at  a 
specific  height. 


j 

EXPINT 


Performs  exponential  interpolations  for 
the  geometry  routines. 


Table  3-4.  Routines  Included  in  Atmospheric  Effects  Module 
(continued) 


Routine 

Description 

SSGEO 

Obtains  attenuator  amounts  from  scatter¬ 
ing  points  along  optical  path  to  the 
extraterrestrial  source. 

PSIDEL 

Calculates  the  relative  azimuth  between 
the  line  of  sight  and  the  direct  solar/ 
lunar  path. 

PSI 

Returns  solar  azimuth  relative  to  line- 
of-sight  at  current  scattering  location. 

DEL 

Returns  solar  zenith  angle  at  any  point 
along  optical  path. 

GEO 

Driver  for  air  mass  subroutines.  Calcu¬ 
lates  attenuator  amounts  for  the  slant 
path. 

SCTANG 

Returns  the  scattering  angle  at  any  point 
along  the  optical  path. 

TRANS 

Calculates  transmittance,  atmospheric 
radiance,  and  solar/lunar  scattered 
radiance  for  slant  path. 

AEREXT 

Interpolates  aerosol  attenuation 
coefficients  to  required  wavenumber. 

HN03 

Determines  nitric  acid  absorption 
coefficient  at  required  wavenumber. 

TRANFN 

Calculates  transmittance  for  ozone, 
uniformly-mixed  gases,  and  water  vapor. 

SOURCE 

Contains  solar  intensity  data  and  calcu¬ 
lates  lunar  intensity. 

TNRAIN 

Calculates  transmittance  of  rain  as  a 
function  of  rain  rate  and  slant  range. 

SSRAD 

Performs  the  layer  by  layer  single 
scattering  radiance  sum. 
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Table*  3-4.  Routines  Included  in  Atmospheric  Effects  Module 
(continued) 


Routine  Description 


PHASEP  Chooses  correct  phase  function  based  on 

relative  humidity,  frequency,  scattering 
angle,  and  model. 

INTERP  Performs  linear  or  logarithmic  interpola¬ 

tion. 

PF  Returns  the  appropriate  phase  function 

from  the  stored  database. 

C1DTA  Returns  water  vapor  band  absorption 

coefficient  at  required  wavenumber. 

C2DTA  Returns  uniformly-mixed  gases  absorption 

coefficient  at  required  wavenumber. 

C3DTA  Returns  ozone  band  absorption  coefficient 

at  required  wavenumber. 

C4DTA  Returns  N2  continuum  absorption 

coefficient  at  required  wavenumber. 

C6DTA  Returns  molecular  scattering  attenuation 

coefficient  at  required  wavenumber. 

C8DTA  Returns  ozone  UV  and  visible  absorption 

coefficient  at  required  wavenumber. 

SLF296  Loads  self-broadened  water  vapor  con¬ 

tinuum  at  296 °K. 


SLF260 


Loads  self-broadened  water  vapor  con¬ 
tinuum  at  260°K. 


Table  3-4.  Routines  Included  in  Atmospheric  Effects  Module 
(continued) 


tinuum . 

MDTA  Model  atmospheric  data. 

PRFDTA  Aerosol  profile  data. 

EXTDTA  Aerosol  extinction  and  absorption  data. 

SF296  Self-broadened  absorption  coefficients 

for  water  vapor  continuum  at  296°K. 

SF260  Self -broadened  absorption  coefficients 

for  water  vapor  continuum  at  260°K. 

BFH20  Foreign-broadened  absorption  coefficients 

for  water  vapor  continuum  at  296°K. 

TRFN  LOWTRAN  transmittance  functions. 

C1D  Water  vapor  band  model  absorption 

coefficients . 

C2D  Uniformly-mixed  gases  band  model  absorp¬ 

tion  coefficients. 

C3D  Ozone  band  model  absorption  coefficients. 

C4D  Nitrogen  continuum  absorption  coeffi¬ 

cients  and  UV  ozone  absorption  coeffi¬ 
cients  . 

PHSDTA  70  averaged  phase  functions  and  truth 

table  identifying  correct  phase  function. 


Routine 


Description 


BATTLE 

AMGREN 

DATRD 
SC  LOS 

BNORM 

MADPR 

UWIND 

RANDOM 

TEMP 

BURN 

MADPUF 

PUFTEM 

BUOY 


Sets  up  geometry,  provides  control  over 
options. 

Calls  routines  for  performing  calcula¬ 
tions. 

Reads  input  deck  from  library. 

Calculates  vectors  for  calculation 
matrix . 

Normalizes  munition  burn  rate  function 
coefficients. 

Calculates  parameters  for  use  by  MADPUF. 

Calculates  average  wind  speed  at  a  given 
height . 

Generates  normally  distributed  random 
numbers. 

Calculates  temperature  at  a  given  height 

Calculates  burn  function  at  a  given  time 

Computes  transport  and  diffusion  of 
random  puffs. 

Calculates  puff  centroid  temperatures. 

Computes  buoyant  rise. 

Calculates  yield  factor  and  acid  mass 
fraction  for  phosphorus  smoke. 


YIELD 


Table  3-5.  Routines  Included  in  Battlefield  Effects  Module 
(continued) 

Routine  Description 

BRISE  Prepares  plume  rise  parameters  for 

exponential  rise  model. 

TRNDF  ACT  II  transport  and  diffusion  model. 

CLDHT  Computes  cloud  temperature  and  buoyant 

rise  in  ACT  II  model. 

ERFUN  Computes  error  function  by  numerical 

approximation. 

SCTCO  Computes  cosine  of  angle  between  two 

vectors. 

BBRAD  Evaluates  blackbody  function. 

SATPR  Computes  water  vapor  partial  pressure. 

CLDRD  Performs  cloud  radiance  calculations. 

PATHR  Computes  radiance  over  a  path  segment. 

SECTR  Computes  sky  sector  radiance. 

Calculates  phase  function  at  a  specific 
angle . 


PHASE 


Table  3-6.  Routines  Included  in  Image  Display  Module 


Routine 

Description 

PRLINE 

Selects  output  choice,  performs  line 
printer  variable  density  plots. 

FIND 

Determines  minimum  and  maximum  values  in 
array. 

COMTAL 

Writes  binary  array  to  a  file,  pads  with 
zeros,  for  reading  by  COMTAL. 

TAPE 

Writes  image  array  to  a  formatted  file. 

Unit  # 

Input  (I) 
Output  (0) 

Description 

Form  of  I/O 

0* 

I 

Battlefield  effects 
module  input 

80  column  card 
images 

3 

0 

Printer  plot  output  for 
producing  variable 
density  printer  plots 

ASCII  characters 
20i644i6  for  37 
gray  levels 

5 

I 

Interactive  (terminal) 
inputs 

List  directed 

6 

0 

Interactive  (terminal) 
outputs  (menus,  prompts) 

Formatted 

7 

0 

Screen  clear  characters 
to  terminal 

Unformatted 

8 

I/O 

Sensor  parameter  library 

Fixed-length 

1233  character 
records 

9 

0 

COMTAL  display  file 

256x256  byte 
data  array 

11 

0 

Formatted  frame  output 

E10.3  format 

15* 

I 

Scene  radiance  map 

80  character 
header.  Each 
data  record 
read  512A4 
format 

oon^ii/ias 


Most  menus  which  allow  the  input  of  values  for  particular 
variables  display  the  current  values  for  the  variables  in 
the  menu  itself.  After  a  new  value  is  input,  the  menu  is 
then  immediately  updated  to  allow  the  user  to  verify  that 
the  value  was  read  in  correctly.  All  numeric  values  are 
read  using  list-directed  format,  which  eliminates  the  neces¬ 
sity  of  entering  the  values  in  a  particular  format.  The 
menus  themselves  accept  single  ASCII  character  inputs.  Only 
the  first  character  entered  is  examined,  and  if  a  character 
is  entered  which  does  not  correspond  to  a  menu  choice,  the 
menu  is  displayed  again  to  allow  the  user  to  make  an  appro¬ 
priate  choice. 

Each  menu  display  is  preceded  by  a  screen  clear  character. 
Since  individual  terminals  require  different  ASCII  charac¬ 
ters  to  clear  the  screen,  the  main  menu  allows  the  user  to 
select  the  appropriate  screen  clear  character  by  choosina 
from  a  library  of  terminal  types.  Currently,  only  screen 
clear  commands  for  the  ADM  31  and  ADM  3A+  terminals  are 

provided.  Appropriate  values  for  the  DEC  terminals  used  at 
TACOM  will  be  added  to  the  library. 

3.1.3.  Order  of  Execution 

The  order  of  execution  of  the  routines  in  the  simulation  may 
be  varied  to  allow  the  user  maximum  flexibility.  However, 
some  care  should  be  taken  to  insure  that  the  individual 
modules  are  executed  in  a  logical  order.  For  example,  it 
would  be  incorrect  to  apply  the  sensor  transfer  functions 

before  running  the  atmospheric  effects  module  to  attenuate  a 
zero-range  radiance  map.  This  section  will  discuss  the 
order  in  which  the  modules  should  be  executed  for  several 

typical  uses  of  the  simulation. 

The  simplest  use  of  the  simulation  involves  reading  in  a 
scene  radiance  map  and  displaying  some  portion  of  it.  To 

accomplish  this,  the  user  should  proceed  as  follows: 

1.  Call  frame  I/O  interface  routines. 

2.  Select  option  to  read  in  a  frame. 

3.  Input  necessary  parameters  for  dimensions  of  the  frame 
to  be  read  in  and  subset  of  the  frame  of  interest. 

4.  Read  in  the  frame  and  copy  a  subset  of  the  frame  to  the 
image  array. 


5.  Return  to  main  menu. 


6.  Select  display  image  option. 

7.  Display  unprocessed  image  in  the  desired  manner. 


To  display  a  scene  radiance  map  attenuated  by  atmospheric 
effects,  the  process  is  only  slightly  more  complicated t 

1.  Follow  steps  1-5  above. 

2.  Call  atmospheric  effects  module. 

3.  Input  appropriate  meteorological  parameters,  path 
geometry  parameters,  and  the  wavelengths  at  which  calcu¬ 
lations  are  to  be  performed. 

4.  Calculate  transmittance  and  radiance. 

5.  Display  image  (steps  6-7  above). 


Instead  of  attenuating  the  scene  radiance  map  with  a  natural 
atmosphere,  attenuation  by  a  battlefield  environment  may  be 
modeled.  The  procedure  is  similar  to  the  previous  case, 
except  that  instead  of  step  4  above,  the  user  should  return 
to  the  main  menu  after  entering  the  atmospheric  effects 
module  parameters  and  call  the  battlefield  effects  module. 
The  user  then  sets  up  the  inputs  for  the  battlefield  effects 
module  and  performs  the  calculations. 

The  effects  of  a  sensor  should  always  be  introduced  after 
the  atmospheric  effects  module  and/or  battlefield  effects 
module  have  been  executed.  The  appropriate  order  of  execu¬ 
tion  is: 

1.  Read  in  frame. 

2.  Call  atmospheric  effects  module  and/or  battlefield  ef¬ 
fects  module. 

3.  Call  sensor  module 

-  input  sensor  characteristics  or  read  in  sensor 
parameters  from  a  library 

select  transfer  functions  to  be  applied 

apply  transfer  functions. 


4. 


Display  image. 


Notice  that  once  a  scene  radiance  map  has  been  read  in,  the 
user  may  process  the  map  in  several  different  ways  without 
having  to  re-read  the  entire  map.  This  is  accomplished 
through  use  of  the  "copy  sub-array"  option  in  the  frame  I/O 
interface  routines.  The  original  scene  radiance  map  is 
always  available  in  memory  and  is  not  modified  by  execution 
of  the  various  modules.  The  user  may  repeatedly  return  and 
copy  some  subset  of  the  radiance  map  into  the  image  array 

for  processing. 

The  user  may  at  any  time  display  the  existing  image  array. 

In  a  given  run,  the  user  may  display  the  original  radiance 

map,  the  radiance  map  attenuated  by  atmosphere,  the  radiance 
map  after  applying  the  battlefield  effects  module,  and  the 
image  after  applying  the  sensor  transfer  functions. 

A  sample  terminal  session  is  provided  in  appendix  B  which 
will  help  the  user  to  become  familiar  with  the  order  of 

execution.  The  session  produces  five  sample  images  which 
include  application  of  the  atmospheric  effects  module, 
battlefield  effects  module,  and  sensor  module. 


3.1.4.  Scene  Geome try 

Scene  geometry  is  specified  by  the  field  of  view  per  pixel 
of  the  scene  radiance  map,  the  angular  dimensions  of  the 
scene  radiance  map,  and  the  range  of  the  scene  radiance  map 
from  the  sensor.  For  running  the  sensor  module  and  the 
atmospheric  effects  module,  no  additional  geometry  infor¬ 
mation  is  required.  For  running  the  battlefield  effects 
module,  the  geometry  is  more  complex  and  additional  input 
parameters  are  required. 


It  is  important  to  note  that  when  running  the  model  to  simu¬ 
late  a  sensor  at  a  range  other  than  the  range  at  which  the 
actual  data  were  taken,  the  field  of  view  per  pixel  must  be 
appropriately  adjusted.  Figure  3-6  demonstrates  the 
geometry.  The  appropriate  relationship  between  field  of  view 
per  pixel  and  range  is: 


F  (r  )  «  (r  /r  J  F  (r  ) 
x  s  ^  m  s  '  x  m 

F  ( r  )  «  (r  /r  )  F  (r  ) 
y  s  ^  m  s  '  y  m 


(9) 


47 


where:  Fx(r),  Fy(r)  »  Field  of  view  per  pixel  at  range 
1  r  from  the  radiance  map  (along- 

scan,  cross-scan) 

rg  =  Range  of  sensor  from  radiance  map  in  the 
simulation 

rm  =  Range  of  sensor  from  radiance  map  when  the 
actual  data  was  taken. 


For  the  battlefield  effects  module,  the  particular  scenario 
of  smoke  munition  locations  and  meteorological  conditions 
must  be  read  in  from  a  library.  Geometry  inputs  which  must 
be  set  up  in  this  library  include  the  coordinates  of  the 
sources  and  a  reference  angle,  PHNOR,  which  specifies  the 
heading  of  north  clockwise  from  the  Y-axis.  The  value  of 
PHNOR  is  used  in  the  model  only  for  the  purpose  of  relating 
wind  direction  to  the  coordinate  system. 

The  values  needed  to  define  the  observer- target  lines  of 
sight  are  input  by  the  observer  in  the  battlefield  effects 
module.  A  grid  of  lines  of  sight  dimensioned  up  to  20x20  is 
specified,  with  the  dimensions  of  the  grid  and  field  of  view 
of  the  grid  defined  by  the  user.  This  grid  is  known  as  the 
smoke  grid.  The  coordinates  of  the  observer  and  the  end  of 
the  lower  left  line  of  sight  in  the  grid  must  be  specified 
using  the  coordinate  system  defined  in  the  input  card  deck. 
This  geometry  is  represented  in  Figure  3-7.  There  must  also 
be  a  way  of  specifying  the  relationship  between  the  smoke 
grid  and  the  scene  radiance  map  array.  To  do  this,  the  user 
must  provide  the  pixel  coordinates  of  the  top  left  corner  of 
the  smoke  grid.  Figure  3-8  describes  the  relationship  be¬ 
tween  the  lines  of  sight  in  the  smoke  grid  and  the  pixels  in 
the  image  array.  Note  that  negative  values  may  be  given  for 
the  pixel  coordinates  of  the  smoke  grid,  although  this  is 
wasteful  of  computer  resources  since  it  means  calculations 
must  be  performed  to  obtain  transmittances  and  radiance 
values  for  lines  of  sight  which  are  outside  the  image 
bounds . 

In  order  to  obtain  the  desired  scene  geometry,  the  user  is 
advised  to  prepare  a  scenario  with  the  smoke  sources  near 
the  origin  of  a  convenient  coordinate  system.  A  qualitative 
drawing  of  the  smoke  source  locations,  wind  direction,  and 
the  resulting  smoke  cloud  should  be  shown  in  a  figure  simi¬ 
lar  to  Figure  3-7.  Then,  the  observer  location  and  smoke 
grid  points  can  be  defined  in  terms  of  this  coordinate  sys¬ 
tem.  To  relate  the  smoke  grid  to  the  image  array,  the  user 
should  display  the  scene  radiance  map  on  the  COMTAL  image 


Wind 

Direction 


Figure  3-7.  Sample  Geometry  Inputs  for  a  Typical  Smoke 

Scenario.  Munition  Locations,  Wind  Direction 
and  PHNOR  are  Specified  in  Smoke  Scenario 
Input  Deck.  Observer  Location,  Specified 
Target  Coordinates,  Dimensions  of  Observer- 
Target  Line  of  Sight  Grid,  and  Field  of  View 
of  Line  of  Sight  Grid  are  Input  via  Battle¬ 
field  Effects  Menu.  In  this  case,  PHNOR  =  40 
Wind  Direction  =  320°. 


processor  before  entering  the  battlefield  effects  module. 
By  noting  the  location  of  the  objects  of  interest  in  the 
image  and  recalling  the  horizontal  and  vertical  fields  of 
view  of  the  image,  the  smoke  grid  may  be  appropriately 
positioned. 

Despite  the  best  efforts  of  the  user,  more  than  one  attempt 
may  be  necessary  to  obtain  the  desired  geometry.  The  simu¬ 
lation  allows  the  user  to  rewind  the  smoke  input  data  file, 
re-copy  the  undistorted  image  from  the  frame  in  memory,  and 
repeat  the  battlefield  effects  module  calculations  using 
different  parameters  without  the  need  to  terminate  the 
program. 


3.1.5  Run-Time  Performance  of  the  Simulation 

The  current  version  of  the  simulation  requires  one  megabyte 
of  virtual  memory  (FFSOO^g  bytes).  This  is  with  the  image 
array  and  the  scene  radiance  map  each  dimensioned  to  256  x 
256.  No  attempt  has  been  made  to  reduce  memory  requirements 
by  operating  the  simulation  in  an  overlaid  mode.  This  will 
not  save  much  memory,  however,  because  more  than  75%  of  the 
memory  in  use  is  required  for  the  scene  radiance  map  and 
image  array. 

Current  CPU  time  requirements  on  the  VAX  11/750  are  a  strong 
function  of  the  paging  demands  on  the  system.  As  a  result, 
time  requirements  do  not  increase  linearly  with  the  number 
of  arithmetic  operations  required.  On  the  University  of 
Michigan's  Amdahl  5860  system,  paging  operations  are  con¬ 
siderably  faster  and  CPU  requirements  vary  nearly  linearly 
with  the  number  of  operations. 

Two  test  cases  were  run  to  demonstrate  the  CPU  time  require¬ 
ments  of  the  simulation.  The  conditions  corresponded  to 
those  listed  in  Table  3-8.  Run  time  and  CPU  time  required 
to  execute  each  of  the  major  routines  are  listed  in  Table  3- 
9  for  operation  on  both  the  VAX  11/750  and  Amdahl  5860  sys¬ 
tems  . 

Of  particular  interest  in  the  time  comparisons  is  the  dra¬ 
matic  increase  in  CPU  time  and  real  time  which  results  when 
the  array  dimensions  are  doubled  and  the  sensor  transfer 
functions  applied  on  the  VAX  system.  Although  the  number  of 
computations  increases  by  only  a  factor  sliqhtly  greater 
than  four,  CPU  time  requirements  increase  by  a  factor  of  13 
and  real  time  requirements  increase  by  a  factor  of  nearly 
25.  These  cases  were  run  with  the  process  limits  described 
in  Table  3-10.  A  total  of  321,587  page  faults  occurred  when 
the  transfer  functions  were  applied  to  the  128  x  128  array 
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Location 

(1,1) 


Figure  3-8.  Relationship  between  Pixel  Array  and  Grid 
Specified  in  Battlefield  Effects  Module. 

For  this  case,  a  9  x  10  Pixel  Array  and  a 
5x5  Smoke  Grid  are  used.  The  Top  Left 
Line  of  Sight  in  the  Smoke  Grid  Corresponds 
to  Pixel  Array  Element  (2,5).  The  Vertical 
Field  of  view  per  Grid  Point  is  Equal  to  the 
Vertical  Field  of  View  per  Pixel;  the  Hori¬ 
zontal  Field  of  View  per  Grid  Point  is  Equal 
to  1.5  Times  the  Horizontal  Field  of  View 
per  Pixel.  Calculations  of  Obscuration  at 
each  Pixel  Array  Location  are  Based  on  Two 
Dimensional  Interpolation  within  the  Smoke 


Table  3-8.  Scenarios  Used  in  Timing  Tests 


SCENARIO  1 

Array  Size:  64x64* 

Sensor  transfer  functions  used:  default 

Atmospheric  Effects  Module  conditions: 

Midlatitude  summer  model  atmosphere 
Rural  haze  model 
23  km  meteorological  range 
0.0  mm/hr.  rain 

800-1240  cm-1  in  20  cm  “1  steps 
Battlefield  Effects  Module  conditions:  none 


SCENARIO  2 

Array  size:  128x128** 

Sensor  transfer  functions  used:  default 

Atmospheric  Effects  Module  conditions: 

Midlatitude  winter  model  atmosphere 

Advection  fog 

1  km  meteorological  range 

1.0  iran/hr.  rain 

800-1250  cm-1  in  5  cm-1  steps 

Battlefield  Effects  Module  conditions: 

4  L8A1  grenades 

5  sec.  after  burst 
Random  puff  mode 


Time  Required  (sec) 


SCENARIO  1 
Initial  Overhead 
Read  Frame 

Sensor  Transfer  Functions 
Atmospheric  Effects  Module 

SCENARIO  2 

Initial  Overhead 
Read  Frame 

Sensor  Transfer  Functions 
Atmospheric  Effects  Module 
Battlefield  Effects  Module 


VAX 

11/750 

Amdahl 

5860 

Real 

CPU 

Real 

CPU 

2 

1.90 

8 

0.276 

12 

11.17 

5 

0.596 

28 

25.52 

4 

1.068 

4 

1.87 

1 

0.105 

2 

1.90 

8 

0.276 

13 

12.49 

5 

0.646 

695 

335.76 

13 

4.491 

6 

4.57 

2 

0.302 

60 

44.36 

18 

2.423 

Priority 


4 


Buffered  I/O  byte  count  quota:  4096 

Timer  queue  entry  quota;  10 

Paging  file  quota:  9851 

Default  page  fault  cluster:  64 

Enqueue  quota:  20 


Working  set  quota: 
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as  compared  with  only  179  for  the  64  x  64  array.  Each 
Fourier  Transform  required  only  five  seconds  of  real  time 
for  a  128  x  128  complex  array,  as  compared  to  an  average  of 
285  seconds  for  a  256  x  256  complex  array. 

3.2.  Sensor  Module  Inputs 

All  inputs  relevant  to  the  sensor  module  are  provided  inter¬ 
actively  by  the  user.  However,  to  free  the  user  from  having 
to  input  all  sensor  parameters  each  time  the  simulation  is 
run,  a  sensor  characteristic  library  feature  has  been  imple¬ 
mented.  Using  this  feature,  the  user  may  recall  any  pre¬ 
viously  defined  sensor  definition  and  then  use  it  or  edit  it 
as  desired.  Similarly,  the  active  sensor  definition  may  be 
saved  in  the  library  for  future  use. 

The  complete  set  of  sensor  menus  is  illustrated  in  Figure  3- 
3.  There  are  several  SENSOR  MENU  options  which  alter  the 
active  sensor  definition.  The  T  option  causes  the  TRANSFER 
FUNCTION  MENU  to  be  displayed.  From  this  menu,  the  transfer 
functions  used  to  model  the  active  sensor  may  be  individual¬ 
ly  selected.  Alternatively,  a  default  selection  of  transfer 
function  "switch"  settings  corresponding  to  any  predefined 
sensor  may  be  loaded  from  the  sensor  characteristic  library. 
The  C  option  in  the  SENSOR  MENU  affords  access  to  the  SENSOR 
CHARACTERISTIC  MENU.  The  L  and  A  options  in  this  menu  con¬ 
trol  access  to  the  sensor  characteristic  library.  A  special 
option  switch  S  controls  the  use  of  transfer  function  switch 
settings  stored  in  the  sensor  library.  If  S  is  false,  the 
transfer  function  switch  table  (TRANSFER  FUNCTION  MENU)  is 
not  loaded  with  the  rest  of  the  sensor  parameters.  Thus, 
TKe"  S  option  will  usually  be  left  in  its  default  true  set¬ 
ting.  The  remaining  menu  options  (  D,W,F,  and  T)  select  the 
menus  from  which  sensor  parameters  are  actually  input.  All 
the  sensor  input  parameters  are  listed  in  Table  3-11;  we 
will  briefly  comment  on  each  menu  in  the  following 
subsections . 

3.2.1  Detector  Characteristic  Menu 

The  inputs  required  by  this  menu  are  relatively  straight 
forward.  Note  that  the  detector  area  is  required  only  for 
noise  calculations,  and  that  the  focal  length  input  units 
are  millimeters.  Also  note  that  options  P  and  S  need  only 
be  selected  if  the  optical  system  blur  and/or  the  LED  trans¬ 
fer  functions  are  used.  All  other  parameters  must  be 
entered  in  all  cases. 
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Table  3-11.  Interactive  Inputs  to  the  Sensor  Module 
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Table  3-11.  Interactive  Inputs  to  the  Sensor  Module  (continued) 
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Table  3-11.  Interactive  Inputs  to  the  Sensor  Module  (continued) 
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SMTFX(I) , 

SMTFY(I)  Arrays  containing  the  along- 

and  cross-scan  stabilization 
MTF  (0-1) 


3.2.2  Spectral  Characteristics  Menu 

There  are  no  optional  inputs  in  this  menu.  A  particularly 
important  parameter  is  the  detector  D*.  This  parameter  de¬ 
pends  on  spectral  wavelength  and  temporal  frequency.  The 
desired  input  is  the  D*  at  the  wavelength  of  peak  spectral 
response,  and  corresponding  to  a  temporal  frequency  that  is 
between  the  1/f  detector  noise  region  and  the  detector  cut¬ 
off  frequency.  Note  that  the  current  version  of  the  sensor 
model  assumes  that  the  detector  noise  is  white  within  the 
modeled  sensor's  electrical  bandpass.  Finally,  note  that 
EWAVE  must  equal  1.0  at  the  wavelength  where  DPEAK  applies. 

3.2.3  Frequency  Characteristic  Menu 

Many  of  the  inputs  allowed  by  this  menu  are  optional.  The 
following  menu  options  must  be  selected  only  if  the  cor¬ 
responding  transfer  function  is  required: 


Transfer  Function 

Detector  -  Temporal 
Electron ics 
Boost 

Along-scan  Sampling 
Stabilization  -  Formula 


Finally,  note  that  the  A  option  allows  limits  to  be  input 
which  are  used  when  the  I  option  (Invert  Transfer  Functions) 
of  the  SENSOR  MENU  is  selected.  The  limits  are  required  to 
prevent  excesssive  noise  amplification  when  performing  aper¬ 
ture  corrections  on  measured  data. 

3.2.4  Tabular  MTF  Inputs 

All  of  the  inputs  requested  by  this  menu  are  optional.  The 
appropriate  inputs  are  required  only  if  the  corresponding 
tabular  transfer  function  is  to  be  used  in  the  sensor 
def in ition . 


3.3  Atmospheric  Effects  Module 


The  atmospheric  effects  module  incorporated  in  the  simula¬ 
tion  is  based  on  the  LQWTRAN6  code.  Major  modifications  to 
the  code  include  the  following: 


1.  Simplification  of  geometry  inputs  to  limit  the  line 
of  sight  to  a  horizontal  (constant-pressure)  path. 
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2.  Replacement  of  all  input  cards  with  a  menu-driven 
calling  routine. 

3.  Elimination  of  the  Navy  maritime  aerosol  module. 

4.  Elimination  of  the  vertical  structure  algorithm. 

5.  Deleting  all  input/output  statements  within  the 
code. 

6.  Elimination  of  the  cirrus  cloud  model. 

7.  Limiting  stratospheric  and  upper  atmospheric  pro¬ 
files  to  default  types. 

8.  Eliminating  the  option  of  the  user  to  input  a  user- 
defined  layered  atmosphere.  The  user  may  still 
input  horizontal  path  parameters  or  a  standard 
model  atmosphere. 

9.  Defaults  are  automatically  provided  for  all  inputs. 


The  following  features  of  L0WTRAN6  were  retained  in  the 
model : 

1.  Scattering  of  solar  radiation  into  the  line  of 
sight. 

2.  Updated  water  vapor  continuum  absorption. 

3.  Relative  humidity-dependent  aerosol  profiles. 

4.  Absorption  and  scattering  due  to  rain. 

5.  Scattering  of  radiation  from  the  ground  into  the 
line  of  sight. 

6.  Calculations  may  be  made  at  any  wavelength  from 
.25  un  to  28  Wi. 


3.3.1.  Methods  of  Operation  of  the  Atmospheric  Effects 
Module 

The  atmospheric  effects  module  treats  natural  obscuration  as 
being  spatially  homogeneous  across  the  scene,  depending  only 
on  range.  In  the  present  version  of  the  simulation,  the 
background  and  target  are  treated  as  being  at  the  same 
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range.  This  will  be  modified  in  the  next  version  to  allow 
separate  background  and  target  images  to  be  combined  at 
different  ranges  from  the  observer. 

There  are  two  different  modes  of  operation  of  the  atmospher¬ 
ic  effects  module  in  the  current  version  of  the  simulation. 
The  first  mode  is  for  the  atmospheric  effects  module  to  be 
called  explicitly  by  the  user  from  the  main  menu.  This 
will,  in  general,  need  to  be  done  at  least  once  whenever  the 
simulation  is  run  in  order  to  set  up  the  meteorological  con¬ 
ditions,  geometry  parameters,  and  wavelengths  at  which  cal¬ 
culations  are  to  be  performed.  After  setting  up  the  inputs, 
the  user  can  either  perform  the  LOWTRAN  calculations  to 
modify  the  scene  radiance  map  or  return  to  the  main  menu  to 
execute  the  battlefield  effects  module. 

The  atmospheric  effects  module  may  also  be  called  implicitly 
from  the  battlefield  effects  module.  In  this  mode,  the  atmo¬ 
spheric  effects  module  is  called  automatically  and  the  user 
has  no  opportunity  to  modify  input  parameters  after  the 
battlefield  effects  module  begins  to  perform  calculations. 
The  battlefield  effects  module  calls  the  atmospheric  effects 
module  twice,  to  compute  transmittance  and  radiance  for  the 
atmospheric  path  between  the  target  and  the  smoke  cloud,  and 
the  path  between  the  smoke  cloud  and  the  observer.  The 
range  used  in  the  atmospheric  effects  module  is  passed 
directly  from  the  battlefield  effects  module  based  on  its 
calculation  of  smoke  cloud  location.  This  range  is  indepen¬ 
dent  of  the  range  which  is  input  by  the  user  in  the  atmo¬ 
spheric  effects  module  path  geometry  menu.  If  the  atmo¬ 
spheric  effects  module  is  called  explicitly  later  in  the 
same  run,  the  user- input  range,  rather  than  the  range  last 
passsed  by  the  battlefield  effects  module,  will  be  used  in 
the  calculations. 

3.3.2  Atmospheric  Effects  Module  Inputs 

The  original  L0WTRAN6  code  required  the  user  to  specify  31 
input  parameters,  with  another  36  input  parameters  required 
if  certain  options  were  selected  by  the  user.  The  modified 
version  of  LOWTRAN 6  included  in  the  atmospheric  effects 
module  requires  only  14  input  variables,  with  three  other 
variables  optional  if  the  user  wishes  to  specify  horizontal 
path  parameters.  Defaults  are  available  for  all  input  vari¬ 
ables.  All  input  variables  may  be  input  from  the  atmospheric 
effects  menu  or  one  of  the  other  menus  called  by  the  atmo¬ 
spheric  effects  menu. 

Table  3-12  describes  the  variables  which  may  be  input  from 
each  of  the  menus  in  the  atmospheric  effects  module.  Every 
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Table  3-12.  Description  of  Options  and  Input  Variables  which  may  be 
Accessed  through  Menus  in  the  Atmospheric  Effects 
Module  (continued) 


Menu 

Name  Option 

L0WTRAN6 

Input 

Variables 

Entered 

Description 

Atmospheric 

Effects  Menu 

M 

Select  meteorological  para¬ 
meter  menu 

P 

Select  path  geometry  menu 

W 

Vl , V2 , DV 

Enter  starting  wavenumber, 
ending  wavenumber,  and  wave- 
number  increment  for  calcu¬ 
lations 

c 

Calculate  transmittance  and 
radiance  using  L0WTRAN6 

0 

Meteorological 
Parameter  Menu 

Return  to  main  menu  without 
performing  calculations 

A 

MODEL* 

Choose  a  model  atmosphere 

H 

IHAZE 

Choose  an  aerosol  model 

V 

VIS 

Enter  visibility  (km) 

R 

RAINRT 

Enter  rain  rate  (mm/hr) 

C 

Enter  refractive  index 
structure  constant  (not 
used  by  L0WTRAN6,  but  used 
by  sensor  module) 

D 

Return  to  default  values 
for  meteorological  inputs 

ti 

Return  to  Atmospheric 

Effects  Menu 

’•JV 


Menu 

Name 

Option 

L0WTRAN6 

Input 

Variables 

Entered 

Description 

Path 

Menu 

Geometry 

H 

HI 

Enter  height  above  ground 

R 

RANGE 

Enter  range  of  path 

J 

I  DAY 

Enter  Julian  day 

T 

TIME 

Enter  time  of  day 

L 

PARMl 

Enter  site  latitude 

A 

PARM2 

Enter  site  longitude 

P 

PS  IPO 

Enter  path  angle  (0=facing 
north,  90  *  east) 

D 

Return  to  default  values 
for  path  geometry  inputs 

Q 

Return  to  Atmospheric 

Effects  Menu 

effort  has  been  made  to  make  it  easy  for  the  user  to  select 
the  appropriate  parameters  for  the  desired  scenario.  The 
appropriate  units  for  the  quantities  to  be  entered  are  in¬ 
cluded  in  the  prompts.  If  further  information  concerning 
inputs  is  required,  the  user  is  advised  to  consult  the  in¬ 
structions  for  using  L0WTRAN6  [1] . 

For  the  user  who  is  already  familiar  with  L0WTRAN6  inputs. 
Table  3-13  describes  which  of  the  L0WTRAN6  inputs  are 
required  by  the  model,  which  are  automatically  set  by  the 
model,  and  which  options  have  been  eliminated.  Default 
values  for  all  inputs  are  provided  in  this  table. 


3.4  Battlefield  Effects  Module  Inputs 


The  battlefield  effects  module  obtains  inputs  both  inter¬ 
actively  from  the  user,  and  from  a  card  image  input  list. 
All  of  the  interactive  inputs  are  found  in  the  Battlefield 
Effects  Menu,  and  concern  the  overall  scenario  geometry. 
The  sole  purpose  of  the  interactive  inputs  is  to  allow  the 
user  to  specify  any  arbitrary  observer-target  viewing 
geometry  easily,  without  re-defining  the  deployment  of 
obscurant  sources.  Information  concerning  specification  of 
the  scenario  geometry  can  be  found  in  section  3.1.4. 


The  card  image  input  list  must  be  prepared  prior  to  running 
the  simulation.  Each  card  image  is  read  using  the  FORTRAN 
format  ( 2A2,6x,7F10 .3) .  The  name  positioned  in  the  first 
columns  of  each  "card"  identifies  the  type  of  information 
contained  in  that  record;  the  card  order  is  generally  unim¬ 
portant.  Table  3-14  lists  the  parameters  which  are  input  on 
each  card.  While  all  these  parameters  are  defined  in 
referenced  reports  [2,3,4],  we  will  comment  briefly  on  each 
card  type  in  the  following  subsections. 


REFD 


The  REFD  card  contains  geometry  and  meteorological  reference 
quantities.  The  variable  PHNOR  is  required  to  define  the 
relation  between  the  scenario  coordinate  system  and  the  com¬ 
pass  headings  used  to  define  wind  direction.  The  parameter 
ZR  defines  the  terrain  roughness  and  thereby  influences  the 
obscurant  vertical  transport.  Pasquill  [11]  discusses  this 
parameter,  however  the  following  table  is  provided  for  easy 
reference : 


Card 

Variable 

Name 

Name 

Code  * 

Default 

Comments 

CARDl 

MODEL 

R 

2  (Mid¬ 
latitude 
summer) 

Model  *  7  not  permitted 

ITYPE 

A 

1 

Horizontal  path  only 

IEMSCT 

A 

2 

Ml 

A 

0 

M2 

A 

0 

M3 

A 

0 

IM 

E 

NOPRT 

E 

T BOUND 

A 

0 

Based  on  model  atmosphere 

SALB 

A 

0.5 

CARD  2 

I  HAZE 

R 

1( Rural) 

Maritime  aerosol  not  per¬ 
mitted 

ISEASN 

A 

0 

Based  on  model  atmosphere 

IVULCN 

A 

0 

ICSTL 

E 

ICIR 

E 

IVSA 

E 

VIS 

R 

23. 

WSS 

E 

WHH 

E 

RAINRT 

R 

0. 

CARD2A 

E 

All  variables  eliminated 

CARD2B 

E 

All  variables  eliminated 

CARD2C1 

E 

All  variables  eliminated 

CARD2D 

E 

All  variables  eliminated 

CARD3 

HI 

R 

0.01 

H2 

E 

ANGLE 

E 

RANGE 

R 

1.0 

BETA 

E 

RO 

A 

0. 

Based  on  model  atmosphere 

LEN 

E 

Table  3-13.  L0WTRAN6  Inputs  and  Defaults  used  in  Atmospheric 
Effects  Module. 


Card 

Name 

Variable 

Name 

Code  * 

Default 

Comments 

CARD3 ' 

HI 

R 

0.01 

P 

0 

1013. 

Only  for  MODEL  -  0 

T 

0 

25. 

Only  for  MODEL  -  0 

DP 

E 

RH 

0 

70. 

Only  for  MODEL  *  0 

WH 

E 

WO 

E 

RANGE 

R 

1.0 

CARD3A1 

IP  ARM 

A 

1 

I  PH 

A 

2 

I  DAY 

R 

270 

I  SOUR 

A 

0 

CARD3A2 

PARMl 

R 

43. 

PARM2 

R 

106. 

PARM3 

E 

PARM4 

E 

TIME 

R 

13.33 

PS  IPO 

R 

90. 

ANGLEM 

E 

G 

E 

CARD3B1 

E 

All  variables  eliminated 

CARD3B2 

E 

All  variables  eliminated 

CARD4 

VI 

R 

800. 

V2 

R 

1250. 

DV 

R 

20. 

CARDS 

E 

All  variables  eliminated 
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Table  3-14.  Card  inputs  to  the  Battlefield  Effects  Module  (continued) 
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1.0  *  Empirical  buoyancy  model 
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done  Card  causes  the  termination  of  battlefield 

module  computations 


i.j  ..  ii  T'W_t  .  .  I 


ZR  (cm) 


Terrain  Type 


0.1 

1.0 

10.0 

100.0 


Open  sand  or  smooth  water 
Short  clipped  grass  in  smooth  terrain 
Open  grass  field  or  agricultural  complex 
Forested  or  urban  environment 


Finally,  the  terrain  scavenging  factor  RC  describes  how  much 
of  the  obscurant  material  which  contacts  the  ground  is  "re¬ 
flected"  back  into  the  air.  The  terrain  is  a  perfect  re¬ 
flector  if  ROl. 

DETD 

The  card  DETD  contains  only  a  single  input,  the  wavelength 
at  which  transmittance  and  radiance  calculations  are  per¬ 
formed  . 

SRCL 

Each  SRCL  card  defines  the  position,  burst  time,  and  muni¬ 
tion  type  for  a  single  obscurant  source.  Up  to  12  separate 
sources  may  be  specified.  Note  that  the  coordinate  system 
is  identical  to  that  used  to  interactively  specify  target 
and  observer  positions  in  the  simulation  itself.  Also,  note 
that  two  munition  type  identification  codes  are  reserved  for 
pre-defined  vehicle  self-screening  grenade  types. 

MET1  and  MET2 

Most  of  the  basic  meteorological  quantities  which  influence 
the  obscurant  transport  and  diffusion  are  input  on  these  two 
cards.  Note  that  either  the  relative  humidity  or  the  dew 
point  temperature  may  be  specified.  The  parameter  WPR  is 
used  to  specify  the  vertical  wind  speed  profile  according  to 
a  power  law  expression  of  the  form: 

U(z)  -  U  (z  -)  (—5 — j  WPR  (10) 

ret  2ref 

where  U  is  the  wind  speed,  z  is  the  height  above  the  ground, 
and  zref  is  a  reference  height. 

Perhaps  the  most  difficult  parameter  to  specify  is  the  in- 
band  surface  irradiance  ESFC.  This  parameter  is  the  net 
irrad lance  from  the  sky  hemisphere  in  a  one  micron  bandpass 
centered  about  WAVE  exclud ing  contributions  from  all  point 
sources  such  as  the  sun  or  moon.  Finally,  note  that  if  the 
ground  is  considered  to  be  a  black  body,  then  the  ground 
reflectivity  ( RFL( 3 ) )  should  be  set  to  0.0. 
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SKYR  and  SKYA 

These  cards  are  used  to  specify  the  relative  radiance  seen 
looking  at  different  points  in  the  sky  hemisphere.  The 
radiance  at  up  to  42  angular  positions  may  be  specified.  Up 
to  seven  radiance  values  may  be  input  on  a  single  SKYR  card; 
three  angular  positions  may  be  input  on  each  SKYA  card.  Of 
course,  the  order  of  angular  and  radiance  inputs  is  critical 
so  that  the  correspondence  between  input  pairs  is  maintain¬ 
ed.  Note  that  the  absolute  magnitude  of  the  SKYR  inputs  is 
of  no  consequence.  The  program  automatically  scales  the 
integral  of  the  diffuse  sky  hemisphere  radiance  to  equal  the 
surface  irradiance  ESFC. 

PPSD 

Up  to  2  PPSD  cards  may  be  used  to  define  point  sources  of 
radiation.  The  parameter  PPSFLX  is  the  surface  irradiance 
in  a  one  micron  bandpass  centered  about  WAVE  due  to  a  point 
source  at  the  specified  angular  position. 

PHFD  and  PHAD 

These  cards  are  used  to  input  the  scattering  phase  function 
of  the  obscurant  media.  Up  to  70  phase  function  values  may 
be  input  on  PHFD  cards  (7  per  card)  while  the  corresponding 
scattering  angles  are  input  on  PHAD  cards.  Only  one  phase 
function  may  be  input;  it  is  used  to  describe  the  scattering 
from  all  obscurant  materials  relevant  to  a  scenario.  There¬ 
fore,  although  any  number  of  different  absorbing  media  may 
be  considered  in  a  single  scenario,  only  a  single  scattering 
media  may  be  correctly  considered. 

SMKD 


This  card  defines  the  relevant  properties  of  the  obscurant 
material.  Note  that  for  the  predefined  smoke  materials 
(phosphorus  smoke  and  EA5763),  only  the  two  optical  para¬ 
meters  applicable  to  the  wavelength  of  interest  (WAVE)  must 
be  input.  The  parameter  OH  quantifies  the  amount  of  heat 
imparted  to  the  obscurant  cloud  in  its  formation;  it  is  used 
to  calculate  the  cloud  temperature  and  buoyancy.  CP  is 
simply  the  specific  heat  of  the  smoke  material.  Finally, 
the  comment  that  only  non-hygroscopic  smokes  may  be  defined 
by  the  user  refers  to  the  fact  that  the  model  does  not  treat 
the  adsorption,  absorption,  and  condensation  reactions  for 
any  material  other  than  phosphorus  smoke,  which  is  pre¬ 
defined  . 


SRCD  and  SGMA 


These  cards  are  used  to  define  smoke  munition  types.  Neither 
card  is  required  if  the  default  grenades  (L8A1  or  XM76)  are 
the  only  source  types  considered.  Note  that  the  product  of 
the  input  munition  fill  mass  and  the  munition  source  effi¬ 
ciency  is  assumed  to  equal  the  total  mass  of  lofted  obscu¬ 
rant  material  produced  by  the  source.  The  munition  source 
efficiency  may  be  set  greater  than  100%  if  the  yield  factor 
of  a  user-defined  hygroscopic  smoke  is  to  be  artificially 
included.  Also,  note  that  the  munition  burn  duration  (RG 
(ITYP,11))  is  described  in  conjunction  with  inputs  provided 
on  the  BURN  card . 

The  inputs  on  the  SGMA  card  refer  to  the  initial  spatial 
distribution  of  smoke  produced  by  the  munition.  A  three 
dimensional  gaussian  mass  distribution  is  used  to  describe 
the  smoke  cloud  produced  by  the  munition  burst,  and  during 
each  successive  modeled  time  step.  The  inputs  are  the  stan¬ 
dard  deviations  of  these  distributions.  A  single  burst  puff 
is  generated  for  each  munition.  Puffs  from  source  material 
which  is  either  scattered  on  the  ground  or  remains  in  the 
munition  canister  are  generated  at  a  rate  of  one  per  second; 
the  smoke  mass  contained  in  each  puff  is  determined  by  the 
munition  burn  function  (see  card  BURN) . 

CNTR 

This  card  specifies  a  number  of  simulation  control  para¬ 
meters.  Assuming  that  STINE  is  specified  as  0.0,  ETIME 
determines  the  time  at  which  the  obscurant  transmission  and 
radiance  map  is  calculated  for  use  in  the  simulation 
scenario.  The  input  DTIME  should  be  set  equal  to  ETIME 
since  the  simulation  ignores  any  intermediate  radiance  map 
calculations . 

The  MODE  input  selects  whether  the  full  single  scatter 
radiative  transfer  solution  is  used.  In  the  fast  radiance 
mode,  the  smoke  cloud  temperature  must  be  input  using  the 
variable  STEMP.  Further,  the  obscurant  cloud  is  treated  as 
a  diffuse  reflector.  When  MODE  is  set  to  1.0,  the  smoke 
temperatures  are  calculated  by  the  model  so  STEMP  need  not 
be  specified.  However,  the  angular  resolution  for  the  sky 
sector  radiance  calculations  must  be  specified.  The  number 
of  sectors  considered  is  the  product  of  NTHE  and  NPHI. 

MAD1 

This  card  allows  either  the  MAD  PUFF  or  the  ACT  II  transport 
and  diffusion  methodology  to  be  selected. 
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TURB 


This  card  allows  parameters  relevant  to  the  MAD  PUFF  trans¬ 
port  and  diffusion  methodology  to  be  input.  If  it  is  not 
known,  the  friction  velocity  may  be  computed  using  the 
following  expression: 

(11) 


USTAR  *  U( z ) *k/( ln( z/ZR) ) 


where  U  is  wind  speed,  z  is  height,  k  is  von  Karman's  number 
(about  0.4),  and  ZR  is  the  terrain  roughness  parameter  input 
on  card  REFD.  This  formula  assumes  that  the  vertical  wind 
profile  is  roughly  logarithmic.  The  remaining  variables 
describe  the  turbulence  velocity  fluctuations  seen  by  an 
observer  moving  with  the  mean  wind.  They  may  be  calculated 
as  follows: 

SIGU  ■  SIG._  *  U 

32  (12) 

SIGV  -  SIGU 

SIGW  =  SIGel  *  U 


where  SIG_Z  and  SIGel  are  the  standard  deviation  of  angular 
wind  fluctuations  measured  in  radians. 

BURN 

The  BURN  card  allows  the  user  to  define  the  rate  at  which  a 
munition  produces  smoke.  The  EO  SAEL  [23]  burn  function 
form  is  used : 

M(t)  *  +  \  B2  +  J  B3  + 

1  Tb  i  i  J  3  Tb  (13) 

1  B4  +  B5  U-0  -  e  86  fc) 

b 


where  Bn  are  the  six  coefficients,  T^  is  the  munition  burn 
duration  input  on  card  SRCD,  and  M(T)  is  proportional  to  the 
cumulative  amount  of  smoke  produced  between  t=0  and  t*T.  For 
convenience,  the  program  normalizes  the  burn  function  so 
that  the  user  need  only  supply  the  correct  relative  weights 
of  the  coefficients;  the  absolute  magnitude  of  the  coeffi¬ 
cients  is  of  no  importance.  Note  the  EO  SAEL  documentation 
[23]  provides  a  convenient  source  of  burn  function  defini¬ 
tions  for  a  wide  variety  of  inventory  and  experimental 
munitions. 
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Finally,  it  should  be  noted  that  a  burn  function  must  be 
input  even  for  the  predefined  munition  types  (L8A1  and 
XM76).  For  convenience,  we  list  burn  functions  for  these 
munitions  below: 


Bum 

Mun ition 

Duration( s) 

B1 

B2 

B3 

B4 

B5 

B6 

L8A1  360.0  0.6  -0.4  0.0  0.0  0.15  0.30 

XM76  1.0  0.0  0.0  0.0  0.0  1.0  999. 


GO  GO  and  DONE 

The  GOGO  card  should  be  placed  following  the  last  data  input 
card;  it  signals  the  program  to  begin  the  smoke  calcula¬ 
tions.  The  last  card  in  the  input  deck  should  be  the  DONE 
card.  This  causes  the  battlefield  module  to  return  control 
to  the  ISSN  after  it  completes  its  calculations.  If  a  DONE 
card  is  not  present,  the  smoke  model  expects  to  see  a  new 
input  deck  for  use  in  an  additional  calculation.  Since  the 
current  version  of  ISSM  only  makes  use  of  the  last  calcu¬ 
lation  performed  by  the  smoke  model,  this  facility  for 
repetitive  calculations  has  no  practical  utility  at  this 
time . 


4.0  Preliminary  Model  Validation  Results 


To  date,  there  has  been  no  attempt  to  validate  the  simu¬ 
lation  as  a  whole,  although  this  is  planned  in  the  future. 
However,  a  preliminary  attempt  has  been  made  to  verify  the 
accuracy  of  the  modules  of  which  the  simulation  is  composed. 
In  the  case  of  the  sensor  module,  we  present  a  comparison 
between  the  measured  and  simulated  performance  of  a  current 
technology  thermal  imaging  sensor.  Although  the  comparison 
is  quite  favorable,  a  considerably  more  thorough  effort  is 
required  in  order  to  truly  validate  the  model. 

In  the  case  of  the  atmospheric  effects  models,  we  have  not 
undertaken  any  effort  to  compare  the  modeled  results  with 
measured  data.  Since  the  natural  and  battlefield  effects 
modules  have  been  developed  and  tested  previously,  we  have 
merely  verified  that  their  performance  in  the  simulation 
matches  the  output  of  the  original  models. 

4 . 1  Sensor  Model  Validation 

In  order  to  provide  a  preliminary  indication  of  the  sensor 
model  validity,  we  have  predicted  the  Minimum  Resolvable 
Temperature  ( MRT )  characteristic  of  an  existing  sensor  for 
which  both  measured  data  and  detailed  modeling  information 
exist.  The  sensor  we  have  modeled  is  a  pre-production  eval¬ 
uation  version  of  the  TADS  thermal  imaging  sensor.  Refer¬ 
ence  24,  which  is  an  unpublished  document  prepared  by  the 
Martin  Marietta  Corp.,  provides  measured  data,  and  NVL 
Static  Performance  Model  inputs  which  describe  the  sensor. 

The  simulation  was  exercised  using  the  referenced  input  data 
to  model  a  sensor  viewing  a  standard  4-bar  MRT  target  pre¬ 
sented  against  a  uniform  background.  We  ran  the  simulation 
for  a  matrix  of  cases  which  included  five  bar  pattern 
spatial  frequencies  and  a  number  of  target  contrast  signa¬ 
ture  levels.  The  resulting  set  of  images  were  then  presented 
to  a  group  of  9  observers.  Within  each  set  of  images  at  a 
particular  spatial  frequency,  the  observers  were  asked  to 
select  the  lowest  contrast  signature  image  in  which  they 
could  both  detect  and  resolve  the  MRT  pattern.  The  observers 
were  allowed  to  optimize  the  displayed  contrast  of  each 
image,  just  as  would  be  the  case  in  a  real  MRT  test.  Also, 
the  observers  were  alloweo  to  adjust  their  viewing  position 
in  any  way  they  desired. 

The  results  of  the  test  are  summarized  in  Table  4-1.  Note 
that  there  is  a  variation  among  the  observers  at  all  spatial 
frequencies  except  for  the  lowest.  In  this  case,  no  observer 
could  detect  the  bar  target  with  a  signature  2.4  times  the 


Table  4-1.  Minimum  Resolvable  Temperature  (MRT)  Simulation  Test  Results 


*.  \  ^  V  V  V  V  v,  v  •- 


V  % 


NVL  predicted  MRT  while  the  target  with  a  signature  4.8 
times  the  NVL  predicted  MRT  was  clearly  visible  to  everyone. 
Interestingly,  the  predictions  based  on  observer  responses 
compare  very  favorably  with  the  measured  data.  Further, 
predictions  made  using  the  NVL  Static  Performance  Model 
directly  seem  to  be  in  relatively  less  agreement  with  the 
measured  data.  However,  while  the  results  of  this  pre¬ 
liminary  test  are  certainly  interesting,  they  cannot  be 
considered  a  sufficient  basis  for  drawing  any  conclusions; 
they  merely  suggest  .that  future  validation  attempts  will 
likely  prove  fruitful. 

4 . 2  Validation  of  the  Atmospheric  Effects  Module 

Validation  of  the  atmospheric  effects  module  was  performed 
by  running  the  module  to  duplicate  the  test  cases  for  hori¬ 
zontal  (constant  pressure)  paths  designated  as  Case  3  in  the 
L0WTRAN5  manual  [10]  and  Case  4  in  the  LOWTRAN6  manual  [1]. 
A  comparison  of  results  is  presented  in  Figures  4-1  and  4-2. 

The  atmospheric  effects  module  was  also  run  to  simulate  the 
conditions  used  to  produce  Figures  40  and  41  in  the  LOWTRAN5 
manual.  These  transmittance  spectra  demonstrate  the  sensi¬ 
tivity  of  the  model  to  variations  in  the  type  of  aerosol 
present.  A  comparison  between  haze  models  is  presented  in 
Figure  4-3,  while  a  comparison  of  transmittance  using  the 
advection  and  radiation  fog  models  is  presented  in  Figure  4- 
5.  The  corresponding  figures  from  the  LOWTRAN5  manual  are 
reproduced  as  Figures  4-4  and  4-6. 

As  can  be  seen,  the  atmospheric  effects  module  duplicates 
LOWTRAN6  results  to  within  the  precision  limits  of  the  cal¬ 
culations.  The  validation  efforts  presented  above  do  not 
validate  the  LOWTRAN6  model,  however,  but  only  demonstrate 
that  the  modifications  to  L0WTRAN6  did  not  affect  the  re¬ 
sults  of  the  calculations. 

Perhaps  the  most  important  deficiency  in  the  LOWTRAN  code  is 
its  inability  to  consider  scatterina  of  radiation  from  all 
sources  into  the  line  of  sight.  In  the  case  of  a  dense 
aerosol,  a  considerable  amount  of  radiation  may  be  scattered 
into  the  line  of  sight  from  the  surrounding  atmosphere.  For 
long  atmospheric  paths  through  a  dense  aerosol,  the  path 
radiance  calculated  by  LOWTRAN 6  may  actually  decrease  with 
increasing  distance.  This,  of  course,  is  physically  incor¬ 
rect.  Table  4-2  demonstrates  this  problem.  This  table  pre¬ 
sents  the  results  of  calculations  for  a  midlatitude  summer 
atmosphere  with  a  radiation  fog  having  a  0.1  km  meteorologi¬ 
cal  range.  Calculations  are  performed  at  800  cm"*  for 
ranges  of  0.05  km  to  10  km.  The  path  radiance  actually  de¬ 
creases  as  the  range  increases  beyond  0.20  km.  As  transmit- 
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Figure  4-2.  Comparison  of  Transmittance  as  Calculated  by  Simulation  and  Test 
Case  4  in  LOWTRAN6  Manual  [1].  Conditions  Correspond  to  a  .3  km 
Path  at  0  km  Altitude,  23  km  Rural  Haze,  10  mm/hr  Rain. 
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Figure  4-3.  Transmittance  Spectra  through  Haze  Models 
Predicted  by  Simulation 
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Figure  4-4.  Transmittance  Spectra  for  a  10-km  Horizontal 
Path  at  Sea  Level  for  the  Rural ,  Maritime, 
Urban,  and  Tropospheric  Aerosol  Models  using 
the  US  Standard  Model  Atmosphere  and  a  Visibil 
ity  of  23  km.  Calculated  by  L0WTRAN5  110] . 
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Figure  4-5.  Transmittance  Spectra  through  Fog  Models 
used  in  Simulation 
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Figure  4-6. 
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Transmittance  Spectra  for  the  Advection  Fog 
(Fog  1)  and  the  Radiation  Fog  (Fog  2)  Models 
for  a  0.2-km  Horizontal  Path  at  Sea  Level, 
with  the  US  Standard  Model  Atmosphere  and  a 
1-km  Visibility,  from  400  to  4000  cm-i.  Cal 
culated  by  L0WTRAN5  [10]. 
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Table  4-2.  Transmittance  and  Path  Radiance  as  a  Function  of 
Range  through  9  Scattering  Medium 


CONDITIONS; 


0.1  km  Visibility 
Radiation  Fog 

Midlatitude  Summer  Model  Atmosphere 

10  m.  Altitude 

800  cm“*  Wavelength 


Range 

Transmittance 

Path 

(W/cm2- 

Radiance 
ster- ym) 

.05 

km 

0.4094 

6.14 

X 

10“6 

.10 

km 

0.1706 

8.43 

X 

10-6 

.15 

km 

0.0714 

9.16 

X 

10"6 

.20 

km 

0.0300 

9.23 

X 

10-6 

.30 

km 

0.0053 

8.76 

X 

10~6 

.40 

km 

0.0009 

8.20 

X 

vo 

1 

O 

H 

.50 

km 

0.0002 

7.72 

X 

10"6 

1.00 

km 

0.0000 

6.58 

X 

10-6 

2.00 

km 

0.0000 

6.24 

X 

10~6 

5.00 

km 

0.0000 

6.22 

X 

10’6 

10.00 

km 

0.0000 

1.24 

X 

10"5 

tance  approaches  the  machine-limited  minimum  value  around 
10“^9,  path  radiance  jumps  to  the  correct  value. 
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As  a  general  rule,  the  atmospheric  effects  module  will  not 
significantly  underpredict  path  radiance  if  transmittance 
due  to  scattering  is  greater  than  10%.  If  extinction  is 
dominated  by  absorption  rather  than  scattering  (as  is  the 
case  with  most  aerosol  models  in  the  8-12  wn  region),  or  if 
the  meteorological  range  is  comparable  to  or  greater  than 
the  range  of  the  path,  the  user  can  be  reasonably  confident 
that  path  radiance  will  not  be  significantly  underpredicted. 

4 . 3  Battlefield  Module  Comparison  Runs 

The  ACTNAD  battlefield  effects  module  was  exercised  in  the 
same  manner  as  LOWTRAN6  in  order  to  verify  that  the  version 
operating  within  the  simulation  provides  the  same  results  as 
known  operating  versions  of  the  model.  Three  cases  were  run 
which  considered  both  the  L8A1  and  XM76  default  munitions. 
Both  of  the  radiance  calculation  options  (set  using  the  MODE 
variable)  were  used.  All  runs  were  made  using  the  MAD  PUFF 
transport  and  diffusion  option. 

The  result  of  the  comparison  was  simply  that  the  simulation 
and  the  original  model  versions  produce  nearly  identical 
results.  In  cases  where  the  random  smoke  transport  had 
little  effect,  the  differences  in  computed  values  were  on 
the  order  of  10”^  relative  to  the  computed  value  itself. 
Note  however  that  since  different  random  number  generators 
are  used  in  the  two  versions,  a  comparison  is  difficult  when 
turbulent  transport  effects  are  significant.  Nevertheless, 
we  believe  the  model  is  functioning  properly.  Of  course,  it 
should  be  emphasized  that  there  is  no  evidence  to  suggest 
that  the  battlefield  obscuration  model  is  valid  for  the 
purpose  of  incorporating  obscurant  clouds  in  observer- 
interpreted  imagery.  This  type  of  validation  must  be  left 
for  future  efforts. 
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